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ABSTRACT 



Context. The inner layers of circumstellar envelopes around AGB stars are sites where a variety of processes such as thermochemical 
equilibrium, shocks induced by the stellar pulsation, and condensation of dust grains determine the chemical composition of the 
material that is expelled into the outer envelope layers and, ultimately, into interstellar space. 

Aims. We aim at studying the abundances, throughout the whole circumstellar envelope of the carbon star IRC +10216, of several 
molecules formed in the inner layers in order to constrain the different processes at work in such regions. 

Methods. Observations towards IRC +10216 of CS, SiO, SiS, NaCl, KC1, A1C1, A1F, and NaCN have been carried out with the IRAM 
30-m telescope in the 80-357.5 GHz frequency range. A large number of rotational transitions covering a wide range of energy levels, 
including highly excited vibrational states, are detected in emission and serve to trace different regions of the envelope. Radiative 
transfer calculations based on the LVG formalism have been performed to derive molecular abundances from the innermost out to 
the outer layers. The excitation calculations include infrared pumping to excited vibrational states and inelastic collisions, for which 
up-to-date rate coefficients for rotational and, in some cases, ro- vibrational transitions are used. 

Results. We find that in the inner layers CS, SiO, and SiS have abundances relative to H 2 of 4 x 10~ 6 , 1.8 x 10~ 7 , and 3 x 10~ 6 , 
respectively, and that CS and SiS have significant lower abundances in the outer envelope, which implies that they actively contribute 
to the formation of dust. Moreover, in the inner layers, the amount of sulfur and silicon in gas phase molecules is only 27 % for S and 
5.6 % for Si, implying that these elements have already condensed onto grains, most likely in the form of MgS and SiC. Metal-bearing 
molecules lock up a relatively small fraction of metals, although our results indicate that NaCl, KC1, A1C1, A1F, and NaCN, despite 
their refractory character, are not significantly depleted in the cold outer layers. In these regions a few percent of the metals Na, K, and 
Al survive in the gas phase, either in atomic or molecular form, and are therefore available to participate in the gas phase chemistry 
in the outer envelope. 
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1. Introduction 

The carbon star envelope IRC +10216 is among the best 
studied astronomical objects (Morris 1975; |Glassgold] 1996; 
|Mauron & Huggins| 1999; ICernicharo etaLl 2000; IMillar"eTaT1 
2000). The central low-mass AGB star CW Leo, losing mass 
at a high rate (1-4 x 10~ 5 M yr _1 ), is embedded in a dense and 
nearly spherical circumstellar envelope composed of molecular 
gas and dust particles. Being a relatively nearby source (1 10-170 
pc), it is one of the brightest infrared objects and richest molecu- 
lar sources - with more than 80 molecules detected - in the sky. 
IRC +10216 is the prototype of carbon star and a challenging 
source for chemical models. 

Molecules observed in IRC +10216 form either in the at- 
mosphere of the central AGB star, from where they are ex- 
pelled (parent species) or in situ in the outer envelope (daugh- 
ter species). Interferometric maps show that molecular emis- 



sion peaks either on the star or farther out in a hollow shell 
of radius 10-20" dLucas et al.l 1995; IGuelin et all 1997). The 
parent species are formed in the dense and hot atmosphere at 
thermochemical equilibrium (TE), although the formation of 
dust grains, shocks, and the penetration of interstellar ultraviolet 
photons through the clumpy envelope ( Agundez & Cernicharo 
2006; Cherchneff 2006; Agundez et al. 2010) can make molecu- 



* Based on observations carried out with the IRAM 30-meter tele- 
scope. IRAM is supported by INSU/CNRS (France), MPG (Germany), 
and IGN (Spain). 



lar abundances to deviate from TE. The chemical composition in 
the inner layers is key as this material, once incorporated in the 
outflowing wind, constitutes the basis for the rich photochem- 
istry and ion-neutral chemistry that takes place in the outer en- 
velope layers. 

Probing molecular abundances in the warm (1000-2000 K) 
and dense (10 8 — 10 14 crrT 3 ) inner layers requires the observation 
of high energy rotational or ro-vibrational lines that lie in the 
submillimeter and infrared regions of the spectrum, which can 
hardly be observed from the ground. Efforts to this end were un- 
dertaken by |Keady & Ridgway| (1993) and |Boyle etaLl ( 1994), 
who observed with ground-based telescopes absorption lines due 
to ro-vibrational transitions of C 2 H 2 , CH 4 , SiH 4 , CS, SiO, NH 3 , 
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and SiS near 10 and 13.5 ^m. lFonfria et al.l (2008) made a sensi- 
tive spectral survey in the 11-14 fim range that unveiled a forest 
of ro-vibrational lines arising from C2H2 and HCN. A recent in- 
terferometric spectral survey in the 0.9 mm atmospheric window 
bv lPatel et al.l (201 1) allowed to detect a good number of narrow 
unidentified lines arising from the innermost envelope, inside the 
acceleration region, that they tentatively assign to molecular ro- 
tational transitions within excited vibrational states. 

Space missions now freely explore spectral regions not ac- 
cessible from the ground. ISO observed a full infrared scan 
from 2.4 to 197 jjm, albeit with a low spatial and spectral res- 
olution, that probed the inner layers of IRC +10216 through 
the lines of abundant molecules such as CO, C2H2, and HCN 
dCernicharo et al.l 1996, 1999). More recently, spectral surveys 
carried out with Herschel with a high angular resolution have de- 
tected many high energy rotational transitions of molecules such 
as CO, HCN, CS, SiS, SiO, and SiC 2 , at low spectr al resolution 
using SPIRE and PACS in the 55-672 fxm range dDecin et al.l 
2010a), and at high sp ectral resolution in the 488-1901 GHz 
range covered by HIFI dCerni charo et al]2010). 

Here we present an exhaustive analysis of the abundance 
and excitation conditions of several molecules formed in the in- 
ner layers of IRC +10216, based on high sensitive observations 
within the full frequency coverage of the IRAM 30-m telescope 
(80-360 GHz). We focus on the parent molecules CS, SiO, SiS, 
NaCl, KC1, A1C1, A1F, and NaCN. The cases of CO, HCN, and 
SiC2, also formed in the inner layers and reachable by the IRAM 
30-m telescope, have been treated recently dDe Beck et al. 1 2012; 
IFonfria et al.l 2008: ICerrncharo et al.1 2010. 2011). The observa- 
tion of a large number of rotational transitions covering a wide 
range of excitation energies coupled to radiative transfer calcu- 
lations allows us to derive accurate abundances throughout the 
envelope. Rotational transitions of CS, SiO, and SiS in highly 
excited vibrational states are observed; for SiS up to v=5, with 
energy levels above 3600 cirT 1 (5200 K). They put severe con- 
straints on the abundance of these species in the innermost re- 
gions. Together with low energy transitions, the abundance can 
be tracked from the inner to the outer layers, allowing to evaluate 
the effects of shocks, grain formation, and gas phase chemistry. 

2. Astronomical observations 

The mm-wave observations of IRC +10216 on which the present 
study is based were made with the IRAM 30-m telescope on Pico 
Veleta (Spain). They cover the main atmospheric windows, at 
A3 mm (80-116 GHz), XI mm (129-183.5), and ^1 mm (197- 
357.5). The observations at 3 and 2 mm were carried out be- 
tween 1986 and 2008 using the old ABCD SIS receivers. A large 
part of the XI mm data have been already presented as a spectral 
survey dCernicharo et al.l 2000). At A3 mm, the most sensitive 
observations were obtained after 2002 in the context of a co- 
herent line survey not yet published (Cernicharo et al. in prepa- 
ration). Part of the spectra at 1 mm were obtained from 1999 
to 2008 using the ABCD receivers although most of the data, in 
particular all spectra with frequencies above 258 GHz, were ob- 
served from December 2009 to December 2010 with the new 
low-noise wide-band EMIR receivers, in the context of a line 
survey (Kahane et al. in preparation). 

Both the old ABCD and the new EMIR receivers are dual 
polarization and operate in single side band. Image side band 
rejections are >20 dB at 3 mm and 2 mm, ~10 dB at 1 mm 
with ABCD receivers, and >10 dB at 1 mm with EMIR receivers. 
Identification of image side band lines during the observations 
with ABCD receivers was done through shifts in the local oscil- 



lator frequency. This was not necessary with EMIR receivers be- 
cause of the markedly different rejection of the horizontal and 
vertical polarizations. The backends provided spectral resolu- 
tions between 1 and 2 MHz. The coarser resolution was used 
above 258 GHz where it corresponds to <2.3 km s , still good 
enough to resolve the lines in IRC +10216's spectra which have 
typical widths of 20-30 km s . 

The observations were made in the wobbler switching ob- 
serving mode by nutating the secondary mirror by ±90" or 
±120" at a rate of 0.5 Hz. The intensity in Figures and Tables 
is expressed in terms of T* A , the antenna temperature corrected 
for atmospheric absorption and for antenna ohmic and spillover 
losses. T* A can be converted into T MB (main beam brightness 
temperature) by dividing by B e ff/F eff , where F e ff, the telescope 
forward efficiency, is 0.95 at 3 mm, 0.93 at 2 mm, 0.91 at 1.3 
mm, and 0.88 at 1 mm for ABCD receivers, and 0.91 at 1.3 mm 
and 0.84 at 0.9 mm for EMIR receivers, and B e ff, the beam effi- 
ciency, is given by B e ff = 0.828 x exp {-(v/341) 2 } for receivers 
ABCD and by B eff = 0.865 x exp {-(v/365) 2 } for EMIR receivers, 
where v is the frequency expressed in GHz. 

The pointing and focus of the telescope were checked every 
1-2 h on nearby planets or the quasar OJ 287. The error in the 
telescope pointing is estimated to be 1-2". The main beam of the 
IRAM 30-m telescope is given by HPBW(") = 2460/v(GHz), so 
that it ranges from 30" at 80 GHz to 7" at 360 GHz. 

IRC +10216 is a variable source with a period of 649 days 
and a variation in the infrared (IR) flux from maximum to min- 
imum phase of 2.2 mag at 1.2 /jm and 0.6 mag at 18 jum 
(ILe Bertrel l992). Since our observations were carried out over a 
time span of several years, we checked for possible time vari- 
ation in the intensity of some lines, especially those that are 
radiatively excited. Most of the studied lines in the A3 and XI 
mm bands were observed at various epochs corresponding to 
different IR phases. As previously discussed in lCernicharo et aLl 
(2000), we have found no evidence of time variability of the line 
intensities that exceeds the calibration uncertainties and could 
be ascribed to changes of the IR flux. Calibration uncertain- 
ties, mostly linked to corrections of atmospheric absorption and 
pointing errors, are estimated to be < 10 % at 3 mm, <20 % at 
2 mm, and <30 % at 1.3 mm. Most of the lines at wavelengths 
shorter than 1.1 mm were observed just once, so that no con- 
straints can be put on their time variability. 

Observations of CS cover the 7 = 2-1 to 7 = 7- 6 rota- 
tional transitions of the v = state, some rotational transitions of 
the excited vibrational states v = 1, 2, 3, as well as several transi- 
tions of the rare isotopomers 13 CS, C 34 S, and C 33 S (see Table[T}. 
For SiO, we observed the 7 = 2-1 through 7 = 8-7 rotational 
transitions of the v = state and 3 transitions of the v = 1 state, 
as well as several rotational transitions of the rare isotopomers 
29 SiO, 30 SiO, Si ls O, and Si n O (see Table^i. In the case of SiS, 
rotational transitions from 7 = 5- 4 to 7 = 19-18 are observed 
in vibrational states up to v = 5. A good number of lines of 
the ground and excited vibrational states were also observed for 
the rare isotopomers 29 SiS, 30 SiS, Si 34 S, and Si 33 S (see Table^. 
Finally, rotational transitions of several metal-bearing molecules 
were observed in their ground vibrational states, including the 
halides NaCl, KC1, A1C1, and A1F, and the cyanide NaCN. For 
this latter, rotational transitions with upper quantum numbers 
7 = 5 to 7 = 23 and K a = to K a = 10 were clearly detected, 
while in the cases of the chlorine-containing molecules both the 
35 C1 and 37 C1 isotopomers were detected (see Tables |4] [5] [6] [7] 
and|8jl. 
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Table 9. IRC +10216's model parameters _ 3 

n(cm ) 2 5 



r(R.) 
10 



20 



50 



T k (K) 



Parameter 


Value 


Distance (D) 


130 pc 


Stellar radius (R,) 


4 x 10 13 cm 


Stellar effective temperature (T,) 


2330 K 


Stellar luminosity (L») 


8750 L Q 


Stellar mass (M„) 


0.8 M Q 


End of static atmosphere (Rg) 


1.2 R, 


Tlnct fntinpncntinn ccirliiit; if? \ 
LJll&i CLHlUCllftdLUJll lclQlLlS f 




End of dust acceleration region (R w ) 


20 R. 


Gas expansion velocity (v cxp ) " 


14.5 km s _1 


Microturbulence velocity (Av tur b) b 


1 km s 


Mass loss rate (M) 


2 x 10- 5 M yr" 1 


Gas kinetic temperature (T k ) 


T, X (r/R,)-°- 55 


Gas-to-dust mass ratio ( — ) 
Dust temperature (7j) 


300 


800 K x (r/R c )-°- 375 



v exp = 5 km s for regions inner to R c , 11 km s 



acceleration region between R c and R w , and 14.5 km s 
b Av lllrb = 5 km s _1 x (R,/r) for regions inner to R c and 1 km s 
the rest of the envelope. 



_1 for the dust 
1 beyond R w , 
in 



,.-0.55 



,.-0.85 



200 R„ and T k oc r 140 beyond 200 R, 



3. Physical model of IRC +10216 

The physical model consists of a spherical envelope of gas 
and dust expanding around a central AGB star. The param- 
eters adopted are given i n Table |9| The assu med distance to 
IRC +10216 is 130 pc dMen'shchikov et all 2001). We con- 
sider an AGB star with an effective temperature T t of 2330 K 
(Rid gway & Keady| 1988), a radius R t of 4 x 10° cm (0.021", 
which is consistent with the value of 0.022" derived from 
infrared interferometry by IMonnier et al.l 2000). The adopted 
values for T* and R* imply a stellar luminosity of 8750 L Q , 
within the range 5200-13,000 L given by Men'shchikov et akl 
(2001). The adopted stellar mass is 0.8 M Q dWinters et al.l 1994; 
IMen'shchikov et al.1 2001). 

The dust condensation radius R c is located at 5 R* and the 
end of the dust acceleration region R w at 20 R„, with the gas 
expansion velocity being 5 km s for the regions inner to R c , 
11 km s _! in the dust acceleration region, and 14.5 km s _1 be- 
yond R w (see Fig. [TJ, in agreement with IR studies of the inner 
layers of IRC +10216 ( [Ready & Ridgwayj 1993; IFonfria et all 
2008). The adopted microturbulence velocity Av tur b is 5 km s _I 
at the stellar photosphere, decreasing down to 1 km s at R c 
(Keady et al. 1988; see Table|9]l. Beyond R c we consider a con- 
stant microturbulence velocity of 1 km s -1 , within the range of 
values 0.65-1.5 km s derived in the literature (Sk inner et al.l 
1999; IDe Beck et al.1 2012). 

The mass loss rate and gas kinetic temperature in the en- 
velope are derived from the modeling of the rotational lines 
/= 1 - to /= 16 - 15 of 12 CO an d 13 CO observe d with the 
IRAM 30-m telescope and HIFI (see lDe Beck et al.l 2012). We 
used a non-local radiative transfer code dDaniel & Cernicharol 
2008), and adopted an abundance of CO relative to H2 of 6 x 
10 4 dCrosas & Mentenl 1997 ; ISkinner et all 19 99) and an abun- 
dance ratio 12 CO/ 13 CO of 45 dCernicharo et al.1 2000). We derive 
a mass loss rate of 2 x 10~ 5 M yr -1 , and a gas kinetic tempera- 
ture that varies with radius as r _a55 for regions inner to 75 R„, as 

-1.40 



between 75 and 200 R„, and as r beyond 200 R* (see 



r -0.85 

Fig.0. 

As concerns dust, we consider spherical grains of amorphous 
carbon with a radius of 0.1 /jm, a mass density of 2 g cirT 3 , and 
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Fig. 1. Particle density n, gas kinetic temperature T^, and expan- 
sion velocity v exp as a function of radius in the inner layers of 
IRC +10216. 



the optical properties of lSuhl (2000). The dust continuum is mod- 
eled using a ray-tracing code that propagates the specific inten- 
sities along a set of impact parameters which are subsequently 
convolved with the telescope beam. The model does not include 
scattering, which becomes important only at short wavelengths, 
below 5 fim. From the modeling of the envelope spectral energy 
distribution (see Fig. [2]) we find a gas-to-dust mass ratio of 300 
and a dust temperature (Tj) radial distribution of the form Tj = 
W0Kx{r/R c )-° 315 . 

The variation of the density of gas particles with radius is 
different depending on the region of the envelope. In the static 
stellar atmosphere, that extends from the photosphere up to a 
radius Rq which we take as 1.2 R„, the density is given by hy- 
drostatic equilibrium 



/ r \ a I GM t (m„) / (R t \ (l - a h 



(i) 



where G and k are the gravitational and Boltzmann constants, 
respectively, a is the exponent in the kinetic temperature law (a 
= 0.55), and (m g ) is the mean mass of gas particles (2.3 amu, 
after considering H2, He, and CO). In the dynamic atmosphere 



10' 



10" - 




X (/im) 

Fig. 2. Spectral energy distribution of IRC +10216 as observed 
by ISO (black line; Cernich aro et al.l 1999), as given by pho- 
tometric measurements (blue crossed boxes; values taken from 
ILe Bertrel (1992) and from the IRAS Point Source Catalog), and 
as calculated by the dust model (red line). 
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shock waves induced by the pulsation of the star extend the cir- 
cumstellar material. For this region, that extends from Rq up to 
R c , we utilize the equation proposed by Cherchneff (1992): 



i UM t {m„)(l -T) I t r Yf- l >\\ 
n = n(R Q )exp\ , f , LL h-l \ (2) 

where y is a dimensionless parameter expressing the ratio of 
the shock strength to the local escape velocity (y = 0.89, see 
ICherchneffl 1992). Finally, beyond the dust condensation radius 
R c , the density at each radius in the envelope is given by the law 
of conservation of mass 



M 



(47rr 2 <m„)v exp ) 



(3) 



which results in a density profile varying as r" 2 , with jumps at 
R c and R w , which reflect the jumps in the adopted expansion 
velocity profile (see Fig.[T]i. 

The treatment of the density profiles is similar to that adopted 

in 



probability, evaluated as (1 - e" Tr )/r v for an expanding spherical 
shell ( Castor 1970), with t v being the optical depth in the radial 
direction within the corresponding shell. The contributions of 
both gas and dust are included in the computation of S v and r v . 
The background radiation field arriving at a given shell is com- 
posed of the cosmic microwave background, the stellar compo- 
nent (geometrically diluted and properly propagated through the 
inner dusty shells), and the dust emission arising from the outer 
shells. 

Apart from the radiation at millimeter wavelengths involved 
in the observed rotational transitions, infrared radiation plays an 
important role for the excitation of some molecules through the 
pumping to excited vibrational states. According to the above 
description, the source of this infrared radiation in each shell is 
on one side local (emission from the warm dust contained in the 
shell that contributes to the local source function) and on the 
other side external (arising from the surrounding shells contain- 
ing warm dust heated by the stellar radiation). 



Agiindez & Cernicharo| (2006), except that we have updated 4.1. Molecular data 



some parameters, such as the stellar mass, and we now do not 
use anymore the phenomenological parameter K„ in Eq. (2 i. The 
density scale is fixed by the mass loss rate through Eq. (3 1, and 
afterwards Eqs. ([TJ and Q determine the density gradient from 
R c down to the photosphere. These two latter equations involve 
several uncertain parameters. For example, the stellar mass is 
not accurately known although it can be confined to the range 
0.6-1 M , based on theoretical arguments and on the late AGB 
stage of IRC +10216 (I Winters et all 1994). This relatively small 
error on the stellar mass translates into an uncertainty of a fac- 
tor of 5 for the density at Rq. At this point it is worth to note 
that densities in the inner layers of IRC +10216 are not par- 
ticularly well constrained by observations and differ by orders 
of magnitude between different modeling studies. For example, 
Willacy & Cherchneff (1998) adopted, in their chemical model 
of the inner wind, a density of 3 X 10 11 cm" 3 at 5 R* rising up 
to 4 x 10 14 cm 3 at 1.2 R t , while ISchoier et al.l (2006a). in their 
model aimed at interpreting SiO observations, adopted a value 
of about 2 x 10 8 cm" 3 at 5 R* rising up to 3 x 10 10 cm" 3 close to 
the photosphere. In our model, densities in the inner layers are 
in between the two above cases, with a value of 1.8 x 10 14 cm" 3 
atR t , 8 x 10 11 cm" 3 at R , and (0.6-1.3) x 10 9 cm" 3 atR c . The 
density of particles is therefore more uncertain in the inner lay- 
ers than in the outer envelope, and this uncertainty translates to 
the molecular abundances derived (see more details in Sec.|4.3|>. 



4. Radiative transfer modeling 

The calculation of the excitation and emergent line profiles 
for the different molecules studied here has been done with a 
multi-shell radiative transfer program based on the LVG (Large 
Velocity Gradient) formalism, which is described in Agundez 
(2009). Briefly, the spherical circumstellar envelope is divided 
into several concentric shells, and statistical equilibrium equa- 
tions are solved in each of them. The radiation field J v at the fre- 
quency v of each transition, needed to solve the statistical equi- 
librium equations, is evaluated in each shell through an escape 
probability formalism as 



j v = (i-P)s v +/3i: 



(4) 



The radiative transfer calculations require as input spectroscopic 
and collision excitation data. Spectroscopic data has been ob- 
tained from a catalog developed by one of us (J. Cernicharo), 
which has been previously used in ICernicharo et al.l (2000) and 
will be published soon. As concerns collision data, ideally one 
would need rate coefficients for inelastic collisions with H2 and 
He involving a large enough number of energy levels and cov- 
ering a wide enough range of temperatures. Since this is not of- 
ten the case, some approximations must be considered. In this 
study we have adopted the following rules for all the studied 
molecules. Inelastic collisions with both H2 and He (with an 



where S v is the local source function, I v s is the specific intensity 
of the background radiation field, and /3 is the so-called escape 



abundance of 0. 17 relative to H2; |Asplund et al. 2009) have been 
considered. When no rate coefficients for collisions with H2 were 
available, as is the case of CS, we scaled those calculated for He 
by multiplying them by the squared ratio of the reduced masses 
of the H2 and He colliding systems. When rate coefficients were 
available for both ortho and para H2, an ortho/para ratio of 3 
was assumed. Wherever we needed rate coefficients involving 
rotational levels higher than those included in the quantum cal- 
culations, we computed them by using the Infinite Order Sudden 
(IOS) approximation (Goldfl am et al.1 1977). No extrapolation in 
temperature has however been made outside the temperature va- 
lidity range of the quantum calculations. Finally, the collision 
excitation rates have been systematically computed from the de- 
excitation rates applying detailed balance. More details are given 
below for individual molecules. 

- CS: we have considered the first 50 rotational levels 
within the vibrational states v = 0, 1, 2, 3. The level en- 
ergies and transition frequencies were calculated from the 
Dunham coefficients provided by Mii ller et al.l (2005), the line 
strengths for rotational transitions were computed from the 
dipole moments for each vibrational state fi v= o = 1.985 D, 
fi v=i = 1.936 D (IWinnewisser & Cookl 1968), gy=2 = 1-914 
D (CDM^J, and /u v=3 = 1.855 D ( |L6pez Pineiro & Tip ping 
1987), while for ro-vibrational transitions the Einstein coeffi- 
cients were obtained from the values calculated by Chandra et al. 
(1995) for the P(l) transition of each vibrational band (see e.g. 
Tipping & Chackerian 1981). Rate coefficients for de-excitation 
through inelastic collisions with H2 and He were taken from 
|Lique et aL (2006), who calculated coefficients for collisions 
with He including the first 31 rotational levels and up to tem- 



See http : / /www . astro . uni -koeln . de/cdms/ 
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peratures of 300 K, and from |Lique & Spielfiedel| (2007), who 
included the first 38 rotational levels within the first 3 vibrational 
states and computed collision coefficients for rotational and ro- 
vibrational transitions at temperatures between 300 and 1500 K. 

- SiO: the first 50 rotational levels of the v = and v = 1 
vibrational states have been included. The Dunham coefficients 
given bv lSanz et al.l (2003). the dip ole moments fi v= p = 3.0982 D 
and yu l= i = 3.1178 D measured by |Raymonda et aL (1970), and 
the Einstein coefficient for the v = 1 — P(l) transition calcu- 
lated by IDrira et al.l (1997) were used to compute the line fre- 
quencies and strengths. The collision rate coefficients calculated 
by |Dayou~ & Balanga (2006) for the first 20 rotational levels and 
for temperatures up to 300 K were adopted. For temperatures 
higher than 300 K and for ro-vibrational transitions we adopted 
the collision rate coefficients used for CS. 

- SiS: we have included the first 70 rotational levels within 
the v = 0, 1, 2, 3, 4, 5 vibrational states. The Dunham coeffi- 
cients and the dipole moment of 1 .735 D given by Muller et al. 
(2007), together with the dipole moments for ro-vibrational tran- 
sitions given by Lopez Pineiro et al. (1987) have been used to 
compute level energies and Einstein coefficients. The rate co- 
efficients of rotational de-excitation through collisions with He 
and ortho/para H2 have been taken from the calculations of 
I Vincent et all (2007) and |Klos & Lique| (2008), which extend up 
to temperatures of 200-300 K. For higher temperatures and for 
ro-vibrational transitions the collision rate coefficients calculated 
bv lTobola et al. 1 (2008) were used. 

- Metal halides: the spectroscopic properties of NaCl, 
KC1, A1C1, and A1F are relatively well known dCaris et al.l 
2002, 2004, IHedderich et al.U 993. IHedderich & Bernathl l992). 
The permanent electric dipole moment is very large for 
Na Cl (9.00117 D; iDeLeeuw et al.l 1970) and KC1 (10.269 
D; van Wachen & Dymanus 1967), and noticeably smaller for 
A1C1 and A1F, 1-2 D and 1.53 D, respectively JLldel 1965). 
For A1C1 we adopted a value of 1.5 D, in agreement with 
the estimation of ILidel (1965) and the value previously used 
by ICernicharo & Guelinl (1987). We have not included vibra- 
tionally excited states since no lines involving such states are 
observed in IRC +10216, and since vibrational line strengths 
for those molecules are poorly known. Recently, Gotoum et al 
(2011,2012) have calculated the rate coefficients for A1F de- 
excitation through collisions with He and para H2 up to tem- 
peratures of 300 and 70 K, respectively. We have also used those 
rate coefficients for A1C1, which has a similar dipole moment 
than A1F. For NaCl and KC1, we have adopted the values used 
for SiS, which extend up to higher temperatures. The lack of ac- 
curate collision rate coefficients at high temperatures for these 
metal halides introduce significant uncertainties in the derived 
abundances. 

- NaCN: we consider rotational levels within the ground vi- 
brational state up to J = 29. The rotational constants used to 
compute the level energies were derived from a fit to the line 
frequencies measured in the spectrum of IRC +10216 and in the 
laboratory by Halfen & Ziurys (201 1^] The dipole moment of 
NaCN is calculated to be 8.85 D (CDMsP). Due to the lack of 
calculated collision data for NaCN we use the following approx- 
imate formula to compute the de-excitation rate coefficients y 
due to collisions with H2: 



where ' and " stand for upper and lower level, respectively, and 
y has units of cm 3 s _1 . Eq. {5} gives a first approximation of col- 
lision rate coefficients for asymmetric rotors, resulting in values 
within one order of magnitude of those calculated for asymmet- 
ric rotors such as SO2 dGreenl l995) and HCOJ (Hammami et al.l 
2007). We, nevertheless, note that the use of this approximate 
formula may introduce significant errors in the excitation and 
abundance of NaCN. 

The lack of collision rate coefficients for A1C1, NaCl, KC1, 
and NaCN may introduce important uncertainties in the derived 
abundances. We may in fact question whether the choice of the 
collision rate coefficients for these molecules is more adequate 
than the simple assumption of local thermodynamic equilibrium 
(LTE). In the case of A1C1, the choice of the collision coeffi- 
cients calculated for A1F does not seem a bad approach based on 
the similar dipole moment and mass of both molecules. In any 
case, the observed lines of A1C1 (see discussion in Sec. |4.3) > arise 
from regions where the level populations are very close to ther- 
malization, and therefore both the use of collision coefficients 
and the LTE approach yield similar results. This is not the case 
for NaCl, KC1, and NaCN (see also discussion in Sec. |4.3| l, for 
which the assumption of LTE results in a very bad agreement 
with the observed lines (the calculated low-/ lines are much too 
weak and the high-/ lines are too strong). We are therefore con- 
fident that our choice of the collision rate coefficients for NaCl, 
KC1, and NaCN yields more reliable abundances than the mere 
assumption of LTE. 



4.2. Modeling strategy 

Once the physical parameters of the envelope, described in 
Sec. [3] and the spectroscopic data and collision rate coefficients, 
detailed in Sec. |4.1| have been established, the only remaining 
parameter needed to compute the emergent line profiles is the 
abundance radial profile for each of the studied molecules. 

As a first trial, we start with a constant abundance from 
the stellar photosphere up to the outer layers where molecules 
are photodissociated by interstellar ultraviolet photons. The 
abundance fall off due to photodissociation is calculated with 
a simple chemical model where the photodissociation rates 
are parameterized as a function of the visual extinction Ay. 
Photodissociation rates for CS and SiO (this latter adopted also 
for SiS) are taken from van Dishoeck et al.l (2006), for NaCl 
(adopted also for KC1 and A1C1) it is taken from Ivan Dishoeckl 
(1998), while for A1F and NaCN we have adopted the educated 
guess 10~ 9 s _I exp(-1.7 x A v ). The visual extinction at each ra- 
dius r is derived by the Nu/Ay ratio, 2.7 x 10 21 cirT 2 , derived 
from our adopted dust parameters and gas-to-dust mass ratio, 
where N H is the column density of total hydrogen nuclei from 
the radius r up to the end of the envelope. The intensity adopted 
for the local interstellar radiation field is half that of the stan- 
dard field of IDrainel (1978): see a discussion on this point in 



logy = -10- 0.2 x (/' - /") - 0.4 x (K' a - K'J) 



(5) 



2 We note that the rotational constants in Table 2 of Half en & Ziurysj 
(201 1) are innacurate, as they result in line frequencies which differ by 
several MHz from the measured ones. 



Agundez & Cernicharo (2006). Once the abundance fall-off due 
to photodissociation is determined, line profiles are calculated 
and compared to the observations. Then, for each molecule, the 
initial abundance is progressively varied until the calculated line 
profiles match the observed ones. 

A constant abundance profile does not always yield a sat- 
isfactory fit to the observed line profiles. More specifically, for 
some molecules it is found that the calculated intensity is over- 
estimated for low-/ lines and underestimated for high-/ lines, 
implying a larger fractional abundance in the warm inner layers 
than in the cold outer regions. Such abundance gradient, likely 
caused by condensation on grains, has been calculated through a 
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simple chemical model which considers that from the condensa- 
tion radius R c gas phase molecules may me adsorbed onto grain 
surfaces, with a rate proportional to a sticking coefficient, the 
thermal velocity of the molecules, the geometrical section of the 
dust grains, and the number of dust grains per unit volume. All 
such parameters are fixed by the physical model described in 
Sec. [3] except for the sticking coefficient which is varied for 
each molecule until the calculated low-/ and high-/ line pro- 
files agree with the observed ones. The radial abundance profiles 
obtained in this way show a pronounced gradient after R c and be- 
come nearly fiat at a radius of about 2 x 10 15 cm. This approach 
does not pretend to precisely simulate the condensation process, 
but at least yields a smooth abundance transition between the in- 
ner and the outer layers of the envelope and gives satisfactory 
results. 

Still, with the approach described above, the calculated 
high-/ lines of some metal-bearing molecules show extra emis- 
sion at the line center, compared with the observed line profiles. 
This implies that the innermost regions, where the expansion ve- 
locity is small, contribute too much in the model to the over- 
all line emission. In such cases we have allowed for a step-like 
abundance decrease in the innermost regions, from the photo- 
sphere up to about 3 R„, where thermochemical equilibrium pre- 
vails (see Agundez & Cernicharo 2006). 

The best fit model, and so the final abundance profile, for 
each molecule is chosen on the basis of the best overall agree- 
ment between the calculated and the observed line profiles. The 
obtained abundance profiles are further discussed for individual 
molecules in Sec. 15.11 

4.3. Results of the radiative transfer modeling 

The observed rotational lines of CS, SiO, and SiS are plotted 
in Figs. [3] [4] and 5] respectively, together with the line profiles 
calculated through the radiative transfer model. In Fig. [T2] we 
plot the derived abundance profiles. For these three molecules 
the observed lines cover a wide range of upper level energies, up 
to 5500 K for CS (level v = 3,/ = 7), 1830 K for SiO (level 
v = 1,/ = 7), and 5440 K for SiS (level v = 5,/ = 19), so 
that they sample different regions of the envelope. Lines from 
the ground vibrational state trace the cool mid and outer layers 
while vibrationally excited lines trace the warm material located 
in the innermost regions (r < 5 R„). The reader may note that 
line widths of vibrationally excited states are significantly nar- 
rower than those of the ground vibrational state, indicating that 
the emission in these lines comes from the innermost regions 
where the gas has not been yet accelerated up to the terminal 
expansion velocity of 14.5 km s . 

As shown in Fig. [3] there is a good overall agreement be- 
tween calculated and observed line profiles for CS, except for 
the v = 1 , / = 3 - 2 line whose observed intensity is substan- 
tially larger than calculated. This line is probably affected by 
a maser amplification, as indicated by the high spectral resolu- 
tion observations made by Highberger et al. (2000). These au- 



thors point towards a mechanism based on infrared pumping to 
the v = 1 state, something that is included in our radiative trans- 
fer treatment, but does not produce a maser amplification in this 
line. A more likely mechanism for this relatively weak maser 
could be the overlap of ro-vibrational CS lines with those of 
abundant molecules such as SiO, a mechanism that explain the 
maser emission observed in oxygen-rich AGB stars in several 
rotational transitions of 29 SiO and 30 SiO, and of SiO in the v = 3 
and v = 4 vibrational states (Gonzalez -Alfonso & Cernicharo] 
1997), and in some SiS rotational lines within the ground vi- 




Fig.3. Rotational lines of CS in IRC +10216 as observed with 
the IRAM 30-m telescope (black histograms) and as calculated 
with the radiative transfer model (magenta lines). Continuous 
lines refer to the best model, using the abundance profile shown 
in Fig. [12] while dashed lines correspond to a model with a con- 
stant abundance of 7 x 10~ 7 from 1 R* up to the photodissocia- 
tion region. Note that this latter model underestimates the inten- 
sity of vibrationally excited lines. 



brational state in IRC +10216 (Fonfria Exposito et al. 2006; see 
below). 

The excitation of CS throughout the envelope is dominated 
by collision processes although the absorption of infrared pho- 
tons emitted by the central star and by the warm dust also plays 
an important role. Radiative pumping at 8 fim populates the vi- 
brationally excited states, enhancing the emission in lines such 
as the v = 3/ = 7- 6 (whose levels lie 5500 K above the 
ground state and which otherwise would not be detectable), but 
also, through downward cascades, affects the population of the 
v = rotational levels. Therefore, a proper treatment of infrared 
pumping warrants a correct estimation of the abundance in the 
inner layers (through the vibrationally excited lines) and allows 
also to accurately determine the abundance in the outer enve- 
lope, which would be overestimated by about 30 % if IR pump- 
ing were neglected. According to our analysis, the abundance of 
CS relative to H2 is 4 x 10~ 6 in the inner layers and drops to 7 x 



12 1. The estimated uncertainty 



10~ 7 in the outer layers (see Fig. 
of the abundance is 50 % in the outer envelope, and a factor of 
2 for the regions inner to R e , where densities are more uncer- 
tain. Fig. [3] shows the improvement in the fit to the line profiles 
brought by the two abundance components model. 

In the case of SiO, we have observed the 7=2—1 through 
/ = 7 - 6 lines of the ground vibrational state as well as 3 lines 
of the v = 1 state. As shown in Fig. [4] there is a good agreement 
between calculated and observed line profiles. A slight differ- 
ence is however found in the v = line shapes (calculated are 
parabolic while observed have a slight double-peak character), 
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Fig. 4. Rotational lines of SiO in IRC +10216 as observed with 
the IRAM 30-m telescope (black histograms) and as calculated 
with the radiative transfer model (red lines). 



which may be indicative of a more extended SiO envelope. A 
better agreement is found, rather than by extending SiO to larger 
radii, by enhancing its excitation in the outer layers. A possi- 
ble source for the required additional excitation could be pro- 
vided by shells with enhanced density, known to be present in 
the outer layers of IRC +10216 ( see Mauron & Huggins 1999; 
ICordiner & Millarl 2009 ; IDe Beck et al.1 2012), or by an infrared 
pumping larger than given by our model. We have nevertheless 
not attempted to improve the agreement for SiO, which overall 
is satisfactory enough, by modifying the physical model, since 
this would affect all the other molecules. The excitation of SiO 
levels in the model is very similar to that found for CS, i.e., is 
controlled by inelastic collisions with H2 while radiative pump- 
ing at 8 /im affects the intensities of the v = rotational lines 
by about 30 % and completely dominates the excitation of the 
vibrationally excited states. The derived SiO abundance, 1.8 x 
10~ 7 relative to H2, remains constant from the innermost layers 



up to the photodissociation region (see Fig. 12 1. The uncertainty 
in the abundance is estimated to be 50 % in the outer layers, 
rising to a factor of 2 in the inner regions, which are sampled 
by just a few vibrationally excited lines and where densities are 
more uncertain. 

Silicon monosulphide (SiS) has been observed in a large 
number of rotational lines covering the v = through v = 5 
vibrational states (see Fig. |5J. The analysis of this molecule is 
complicated by the fact that some rotational lines within the 
ground vibrational state (concretely the J = 1 1 - 10, J = 14-13, 



and J — 15-14) show a profile with enhanced emission due 
to maser amplification at certain velocities, which are probably 
caused by the overlap at 13.5 /mi of ro-vibrational transitions 
of SiS with those of C2H2 and HCN (Fonfria Exposito et al. 
2006). The agreement between calculated and observed fine 
profiles is somewhat poorer than in the cases of CS and SiO. 
On the one hand our radiative transfer calculations do not in- 
clude infrared overlaps, something that clearly affects the three 
maser lines previously mentioned but can also have an influ- 
ence on the excitation of other vibrationally excited lines (see 
Fonfria Exposi to et al.| 2006), and on the other hand the observed 
intensity of some high frequency lines may have non negligible 
errors due to calibration (estimated to be around 30 % in the AQ.9 
mm band), to pointing inaccuracies (which are typically <2", al- 
though the main beam of the IRAM 30-m telescope is as small as 
7" at 360 GHz and vibrationally excited lines have an emission 
size of less than 1"), or even due to time variability (difficult to 
evaluate due to the lack of observations at different times in the 
A0.9 mm band). 

Infrared pumping affects little the level populations within 
the ground vibrational state of SiS, but controls, at practically 
the same extent than collisions, the excitation of the vibra- 
tionally excited lines. As in the case of CS, we need to increase 
the abundance of SiS in the inner layers with respect to that in 
the outer envelope. We derive for SiS an abundance relative to 
H 2 of 3 x 10~ 6 in the inner layers, decreasing down to a value 
of 1.3 x 10~ 6 in the mid and outer layers (see Fig. 12 1. The 
estimated uncertainty in the abundance of SiS is a factor of 2 in 
the inner layers and 60 % for the outer envelope. 

The metal halides NaCl, KC1, A1C1, and A1F are all observed 
through rotational transitions within the ground vibrational state. 
No vibrationally excited lines are detected, probably because of 
the low abundance of these species. The observed and calculated 
line profiles of Na 35 Cl and K 35 C1 are shown in Fig. [6] and Fig. [7] 
respectively. These two species have similar dipole moments, 
so that the larger intensity of the NaCl lines largely reflects its 
larger abundance with respect to KC1. The abundance derived 
for Na 35 Cl in the inner layers is 1.3 x 10~ 9 relative to H2 (with 
an uncertainty of a factor of 2), although we find necessary to 
decrease the abundance in the thermochemical equilibrium re- 
gion (1-3 R„) down to a value of 10~ n (see Fig. |l3|> to avoid a 
too stronge mission for the high-/ lines near the central velocity 
(see Fig.|6]l. In the case of K 35 C1 the derived abundance profile 
is somewhat more elaborated (see Fig. 



13 1. KC1 abundance is 



found to have a slight gradient going from 5 X 10 at 5 R* 
down to 3.7 x 10~ 10 in the outer layers, with an estimated error 
of a factor of 2. As with NaCl, the abundance is decreased in the 
1-3 R* region, although in this case by just a factor of 2.5. 

The observed and calculated line profiles of A1 35 C1 and A1F 
are shown in Fig. [8] and Fig. [9] respectively. A1C1 and A1F have 
a much smaller dipole moment than NaCl, but show similar 
line intensities, which suggest they are fairly more abundant. 
Radiative transfer calculations confirm this point and yield abun- 
dances relative to H 2 of 5 x 10~ 8 for A1 35 C1 and 10~ 8 for A1F 
(with an uncertainty of 50 % increasing up to a factor of 2 in the 
regions inner to R c ). In the thermochemical equilibrium region 
(1-3 R„) the abundance of A1 35 C1 is decreased, as for Na 35 Cl 
and K 35 C1, by a factor of 5 (see Fig. 13 I. It is worth to note that 



the observed lines of A1F are slightly more U-shaped than those 
of A1C1. Such behavior can be accounted for to some extent by 
assuming that A1F is more extended than A1C1 due to a lower 
photodissociation rate, as assumed in our model (see Sec. 4.2 1. 
Still the observed line profiles of A1F have a slightly more pro- 
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Fig. 5. Rotational lines of SiS in IRC +10216 as observed with the IRAM 30-m telescope (black histograms) and as calculated with 
the radiative transfer model (blue lines). 



nounced U-shape than calculated (see Fig. [9]). Here, as occurs 
with SiO (see above), a better agreement cannot be obtained by 
extending A1F to larger radii, but by enhancing the excitation 
of A1F in the outer layers. In fact, assuming LTE excitation for 



A1F brings the calculated line profiles in close agreement with 
the observed ones. This is somewhat paradoxical since A1F is 
the only metal-bearing molecule for which collision rate coef- 
ficients have been calculated. This slight discrepancy is, never- 
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Fig. 6. Rotational lines of Na 35 Cl in IRC +10216 as observed 
with the IRAM 30-m telescope (black histograms) and as calcu- 
lated with the radiative transfer model (green lines). Continuous 
lines refer to the best model, using the abundance profile shown 
in Fig. 13 while dashed lines result from a model without the 



abundance decrease in the 1-3 R* region. Note that dashed lines 
overestimate the intensity at the center of high-/ lines. 



theless, unlikely to have a significant impact on the A1F derived 
abundance. 

The excitation of the metal halides in IRC +10216 is likely to 
be dominated by inelastic collisions. The role of infrared pump- 
ing to vibrationally excited levels, not included in our radiative 
transfer calculations, is difficult to evaluate mostly due to the 
lack of ro-vibrational line strengths. To get a rough idea of the 
importance of infrared pumping we have run radiative transfer 
models including the first vibrationally excited state. As dipole 
moment for the v=l— >0 band we have adopted that of SiS (0.13 
D), scaled by the ratio of the halide to SiS permanent electric 
dipole moments. These calculations indicate that infrared pump- 
ing does not significantly affect the intensities of the observed 
lines, except for the high-/ lines of NaCl which are slightly 
more intense with respect to the case when infrared pumping 
is not included. We may conclude that it is unlikely that infrared 
pumping plays an important role in the excitation of the v = 
rotational levels of these metal halides, which is very likely dom- 
inated by collision processes. 
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Fig. 7. Rotational lines of K 35 C1 in IRC + 102 1 6 as observed with 
the IRAM 30-m telescope (black histograms) and as calculated 
with the radiative transfer model (red lines). 



Even in the absence of IR pumping, the ground vibrational 
state populations of metal halides depart from LTE, according 
to our statistical equilibrium calculations. Fig. 10 shows, for se- 
lected rotational transitions, the calculated ratio of the excita- 
tion to kinetic temperature as a function of radius. We see that 
for the 4 metal halides this ratio is equal to 1 in the inner lay- 
ers, i.e. the rotational levels are thermalized, whereas at larger 
radii the decrease in density makes rotational levels to be sub- 
thermally populated yielding excitation temperatures lower than 
the kinetic temperature. NaCl and KC1, with a large dipole mo- 
ment, start to deviate from thermalization at radii much shorter 
than A1C1 and A1F, which due to their moderately low dipole 
moments keep close to thermalization throughout a large part of 
the envelope. This suggests that the choice of the collision rate 
coefficients should not be critical for A1C1 and A1F. 

Finally, as concerns the cyanide NaCN, a fairly large num- 
ber of rotational lines within the ground vibrational state has 
been observed (see Fig. [TT] ). Although the agreement between 
calculated and observed line profiles is overall satisfactory, tak- 
ing into account the large uncertainty of the collision rate coef- 
ficients, the calculated high-/ high-^T fl lines are too strong and 
the low-/ and low— K a lines too weak. We note that the discrep- 
ancy is larger if we replace statistical equilibrium calculations by 
LTE. The adopted photodissociation rate of NaCN is also very 
uncertain (see Sec. |4.2[). Our model points toward a lower pho- 
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Fig. 8. Rotational lines of A1 35 C1 in IRC +10216 as observed 
with the IRAM 30-m telescope (black histograms) and as calcu- 
lated with the radiative transfer model (blue lines). 



todissociation rate than adopted, which would extend the NaCN 
envelope increasing the intensity of low-/ lines, although this 
conclusion is hampered by the uncertain collision rates adopted, 
which are likely to cause most of the observed discrepancies. 
The role of infrared pumping to vibrationally excited states is 
difficult to evaluate due to the lack of spectroscopic and line 
strengths data for such states. The derived abundance for NaCN 
is 3 x 10" 9 relative to H2, with an estimated uncertainty of a fac- 
tor of 3. Decreasing the NaCN abundance in the 1-3 R* region, 
down to a very low value (1CT 12 ), improves somewhat the fit to 
the observed high-/ high-/T fl lines. 

5. Discussion 

5.1. Molecular abundances 

A compilation of molecular abundances for the inner layers of 
IRC +10216, resulting from this work and from previous stud- 
is given in Table 10 The abundance profiles derived in 
T2| for CS, SiO, and SiS, and in 
, A1F, and NaCN. For chlorine- 



les 



this study are shown in Fig 
Fig. [13] for NaCl, KC1, A1C1 
containing molecules, the abundances given in Table 10 and 
shown in Fig. 13 include both the 35 C1 and 37 C1 isotopic species, 
with a 35 C1/ 37 C1 ratio of 2.9, as derived in Sec.^ 
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Fig. 9. Rotational lines of A1F in IRC +10216 as observed with 
the IRAM 30-m telescope (black histograms) and as calculated 
with the radiative transfer model (magenta lines). 



In Fig. 12 and Fig. 13 we also show the abundances com- 
puted under thermochemical equilibrium, which is expected to 
be valid from the stellar photosphere up to a radius of ~3 R*, 
from which the decrease in density and temperature causes an 
increase in the chemical time scale so that the abundances ex- 
perience a quenching effect (see |Agundez & Cernicharo"] 2006). 
The thermochemical equilibrium calculations use solar elemen- 
tal abundances ( Asplund et al. 2009), except for carbon whose 



abundance is increased over the oxygen value by a factor of 
1.5. More details on the thermochemical equilibrium calcula- 
tions can be found in Agundez (2009). By looking at Fig. 12 
and Fig. 13 we note that for most of the studied molecules - a 
notable exception is CS - the abundance calculated under ther- 
mochemical equilibrium shows a positive gradient in the 1-3 R* 
region, i.e. the abundance close to the photosphere is much lower 
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Fig. 10. Calculated ratio of excitation temperature to kinetic tem- 
perature (/ ex / Tk) as a function of radius for selected rotational 
transitions of the metal halides Na 35 Cl, K 35 C1, A1 35 C1, and A1F 
in IRC +10216. 
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Fig. 12. Abundances of CS, SiO, and SiS in IRC +10216, as de- 
rived from the radiative transfer model (continuous lines) and as 
calculated through thermochemical equilibrium in the innermost 
regions of the envelope (dashed lines). A vertical dashed line in- 
dicates the outer boundary where thermochemical equilibrium is 
valid (~3 R«). 
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Fig. 13. Abundances of NaCl, KC1, A1C1 (including 35 C1 and 
37 C1), A1F, and NaCN in IRC +10216, as derived from the radia- 
tive transfer model (continuous lines) and as calculated through 
thermochemical equilibrium in the innermost regions of the en- 
velope (dashed lines). A vertical dashed line indicates the outer 
boundary where thermochemical equilibrium is valid (~3 R*). 



than in the surroundings of the quenching region (~3 R„). This 
behavior justifies the decrease in abundance in the 1-3 R* region 
adopted for various metal-bearing molecules. 

For some species (CS, SiS, and KC1) we have found 
necessary to consider a negative gradient in the abundance, 
with a higher value in the inner layers than in the mid and 
outer regions, in order to adequately reproduce the wide range 
of low- and high-energy observed transitions. The scientific 
rationale for such a gradient is that the studied molecules contain 
refractory elements and are likely to deplete from the gas phase 
to condense onto dust grains. Therefore, we consider that the 
gradient starts at the condensation radius R c and continues 
up to a radius of a few 10 15 cm. Farther out, at about 10 16 
cm all molecules start to be photodissociated by the ambient 
interstellar ultraviolet field. The observation of a large number 
of rotational lines for the different molecules studied, involving 
levels with a wide range of energies (including highly excited 
vibrational levels) allows us to trace the material, and therefore 
to determine the abundance, from the innermost layers out to 
the outer envelope. 

Carbon monosulphide (CS) is found to have an abundance 
relative to H2 of 4 X 10~ 6 in the inner layers, decreasing in the 
mid envelope down to a value of 7 x 10~ 7 at a radius of 2 x 
10 15 cm, which serves as input value for the photochemistry 
taking place in the outer envelope (see Fig. (12) . Our estimate 
for the inner layers may be compared with the results of the 
thermochemical equilibrium calculations, which yield a steep 
abundance gradient around the quenching region, from 2 x 10~ 5 
at 2 R, down to 3 x 10~ n at 4 R*. The intersection of our derived 
abundance profile with that calculated under thermochemical 
equilibrium indicates that the abundance of CS is quenched at 
a radius of 2.4 R„. The model built by Willacy & Cherchneff 
(1998), which considers the effect of shocks induced by the stel- 
lar pulsation on the chemistry, predicts abundances somewhat 
lower in the inner layers, 2-6 x 10~ 7 in the 2.7-5 R* region. 
Our findings for the inner layers are in perfect agreement with 
the previous results of Keady & Ridgway (1993), who derive 
an abundance of 4 X 10~° for the region inner to 17 R* from 



Highberger et al. 
J (2009)" >9.3 x 



observations of ro-vibrational lines in the infrared, but differ 
significantly from the estimates made by 
(2000), 3-7 x 10~ 5 at 14 R„, and bv lPatel eta 
10~ 6 within a radius of 7 R„, which are based on observations 
of rotational lines in vibrationally excited states and an analysis 
involving column densities and mean excitation temperatures. 
As concerns the abundance in the mid and outer envelope, 
our estimate, 7 x 10~ 7 , is somewhat higher than the previous 
value derived by Hen kel etaTI (1985). 1.2 x 10~ 7 , based on 
the observation of low-/ rotational lines within the ground 
vibrational state. 

For silicon monoxide (SiO) we derive an abundance rela- 
tive to H2 of 1.8 x 10~ 7 , from the innermost circumstellar lay- 
ers up to the photodissociation region. Yet, as in the case of 



CS, chemical equilibrium calculations (see Fig. 12 1 yield a steep 



gradient near the quenching region, with a value as low as 3 x 
10~ 10 at 2 R„ increasing up to 10~ 6 at 4 R*. The quenching ra- 
dius for SiO would be located around 3.6 R„, somewhat farther 
than that of CS. We find no evidence of such a low SiO abun- 
dance in the innermost layers, perhaps because only a few vi- 
brationally excited lines tracing this region enter in our analysis. 
Calculations based on chemical kinetics (see model MH in Fig. 3 
of Agundez & Cernicharo 2006) indicate that the abundance of 
SiO increases when moving away from the photosphere due to 
the reaction between Si and CO, which after a few stellar radii 
becomes too slow compared to the dynamical time scale, so that 
the abundance of SiO gets quenched. The chemical model with 
shocks of Willacy & Cherchneff (1998) predicts an abundance 
of 2-7 x 10'' in the 1.9-5 R* region, which is in good agree- 
ment with our derived value. 

As concerns previous observational studies, one of the most 
complete dealing with the abundance of SiO in IRC +10216 is 
that made by Schoi er et al.l (2006a). who used observations of 
several rotational lines within the ground vibrational state, in- 
terferometric observations of the v = 0/ = 5— 4 line, as well 
as ro-vibrational lines observed at 8 fim to derive the abundance 
of SiO throughout the envelope. They found an abundance of 
3 x 10~ 8 from the photosphere up to 1.7 x 10 14 cm, a high 
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abundance component of 1.5 x 10~ 6 extending up to 4.5 x 10 14 
cm (which is necessary to reproduce the ro-vibrational lines ob- 
served), and an abundance of 1.7 x 10~ 7 from this latter radius 
up to the outer layers. The main difference with our results is that 
Scho ier et al.l (2006a) find necessary to adopt a high abundance 
component in the inner layers. A higher abundance than derived 
by us, 8 x 10~ 7 , was also found by Keady & Ridgway (1993) in 
their analysis of the observations at 8 jum. These authors adopt 
lower densities than us for the inner layers, they assume that 
Eq. (|3]l applies down to the photosphere, and this could explain 
the higher abundances they get. On the other hand, the abun- 
dance derived by Keady & Ridgway (1993) for CS is in very 
good agreement with our value (see above), in spite of the very 
different densities adopted for the inner regions. Given the dif- 
ferent nature of the lines used to trace the abundance in the inner 
layers, infrared ro-vibrational lines versus millimeter rotational 
lines in vibrationally excited states, it is difficult to explain why 



the abundances found by |Keady & Ridgway (1993) and by us 
are different for SiO but not for CS. A forthcoming study of 
new high spectral resolution observations (R~ 100,000) of IRC 
+ 10216 around 8 //m obtained from ground (Fonfria et al. in 
preparation) will help to solve the discrepancy found between 



our derived SiO abundance and that found by |Keady & Ridgway 
(1993). 

More recently, De cin et al.l (2010a) have observed a large 
number of SiO rotational transitions in IRC +10216, some of 
them in the v = 1 state, at low spectral resolution with the 
SPIRE and PACS instruments on board Herschel. These authors 
derive an abundance of 10~ 7 for SiO from the innermost re- 
gions up to the outer layers, in good agreement with our findings. 

According to our best model, silicon monosulphide (SiS) is 
present in the inner layers with an abundance relative to H2 of 3 
X 10 6 , decreasing in the mid envelope down to a value of 1.3 x 

the thermo- 
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10" 6 at a radius of 2 x 10 15 cm. As shown in Fig. 
chemical equilibrium abundance of SiS shows a steep gradient, 
rising from 10~ 7 at the photosphere up to a maximum value of 
2.6 x 10~ 5 beyond 2.5 R*. Our derived abundance is substan- 
tially lower than this latter value, which points towards the SiS 
abundance being quenched in a region inner to 2.5 R„, concretely 
at 1.9 R„, as deduced from the intersection of the continuous and 
dashed lines in Fig. [72] Calculations based on chemical kinetics 
i fWillacy & Cherchneff| 1998; |Agundez & Cernicharo| 2006) re- 
sult in SiS abundances of 2-3 x lO -3 in the regions inner to 5 
R„, about one order of magnitude larger than the value we find. 
These models are, however, affected by uncertainties due to the 
poor knowledge of rate constants for reactions involving Si- and 
S-bearing species. 

As concerns previous observational studies, there is a 
general consensus on the presence of a gradient in the abun- 
dance of SiS, with reported values in the range 10 6 — 10 -5 
in the inner layers (>6.5 x 10 " 7 by ITurnerl 1987; 7.5 x 10~ 6 
by |Bieging & Nguyen-Q-Rieu| 1989; 4.3 x \Q 5 close to the 
stellar surface decreasing down to 4.3 x 10~ 6 at 12 R* by 
Boyl e et aT] 1994) and reported values in the range 10~ 7 — 10 6 
in the o uter envelope (2.4 x 10~ 7 by ISahai et all 1984; 5 x 
10~ 7 by |Nguyen-Q-Rieu et al.| 1984; 5 x 10" 7 by IHenkel etall 
1983; L5-4 x 1Q-' by |Henkel et al.l 1985; 6. 5 x 10' 7 by 
Biegi ng & Nguyen-Q-Rieu] 1989; 14 x 10" 6 b y ISchoier et all 
2007). The more recent study of IDecin et al.l (2010a). based 
on observations with SPIRE and PACS that are particularly 
sensitive to the warm inner envelope, reports a SiS abundance of 
4 x 10~ 6 , in good agreement with our finding for the inner layers. 



Table 10. Molecular abundances in the inner layers of IRC 
+10216 



Molecule Abundance relative to H2 Reference 

1-5 R, -> 2 x 10' a cm 

~CS 4 x 1CT 6 ^ 7 x 10- ' (TT 

SiO 1.8 x 10~ 7 (1) 

SiS 3xl0~ 6 -> 1.3 x 10- 6 (1) 

Nad" 1.8 x 10-" (1) 

KC1" 7 x 10-'° -> 5 x 10-'° (1) 

A1C1" 7 x 10~ 8 (1) 

A1F 1 x 10~ 8 (1) 

NaCN"' 3 x IP' 9 (1)_ 

~CO 6 x 10~ 4 (TT 

C 2 H 2 8 x 10~ 5 (2) 

HCN 2 x 10~ 5 (2) 

CH 4 3.5 x 10~ 6 (3) 

NH 3 2 x 10- 6 (4) 

SiH 4 2.2 x 10- 7 (3) 

SiC 2 2 x 10~ 7 (5) 

H 2 1 x 10~ 7 (6) 

HC1 1 x 10- 7 (7) 

HCP 2.5 x 10~ 8 (8) 

C 2 H 4 2 x 10~ 8 (9) 

HF 8 x 10- 9 (7) 

PH 3 8xlO- 9 (10) 

H 2 S 4 x IP' 9 (11) 

Notes: " Abundance in the 1-3 R» region is reduced to 1 x 10~" 
(NaCl), 2 x 10-'° (KC1), 1 x 10~ 8 (A1C1), and 1 x 10~ 12 (NaCN). 
References: (1) this study (NaCl, KC1, and A1C1 abundances in- 
clu de both 35 C1 and 37 C1 isotopomer s); (2) IFonfrfa et al.l (2008); 
(3) |Keady & Ridgway| (19 93); (4) |Hasegawa et ah] (2006); (5) 
ICernicharo et al.l QOlO ); (6)IDecin et al.l (2010b); ( 7) | Agundez etaT 
(2011); (8) based on | Agundez et al | (2007); (9) IGoldhaber et al. 
(1987); (10) based on Agu ndez et alT(2008) ; (11) derived from the 
line 1 10 - 1q,i observed by Cernicharo et al. (2000). 



The most abundant sulfur-bearing molecules in IRC 
+ 10216 are therefore CS and SiS, H2S being a minor species 
(see Table 10 1. Assuming a solar elemental abundance, we 
conclude that 27 % of the total available sulfur is locked in gas 
phase CS and SiS molecules in the inner layers. The fraction 
decreases down to 7.6 % in the outer envelope, an evidence of 
sulfur depletion onto dust grains in the dust formation region. 
A fraction of the remaining sulfur could be in atomic form, 
but it is likely that the bulk of sulfur is in the form of solid 
condensates such as MgS ( |Goebel & Moseleyj 1985). Several 
silicon-bearing molecules such as SiS, SiO, SiC2, and SiH 4 
are present in the inner envelope, all them locking up to 5.6 % 
of the total available silicon. This small fraction, together with 
the fact that SiS (the most abundant Si-containing molecule) 
gets depleted by a factor of 2.3 after the dust condensation 
region, indicates that silicon massively condenses into dust 
grains, mainly forming SiC ( ITreffers & Cohenl 1974). We see 
no evidence of SiO depletion at the dust formation region in our 
SiO abundance profile. One would expect a depletion similar to 
that of SiS, silicon being a refractory element that is thought to 
easily condense onto dust grains in the envelopes of AGB stars 
(ISchoier et al.l 2006b). This may result from our small sample 
of vibrationally excited SiO lines, which limits our sensitivity to 
the innermost layers. 

The metal halides NaCl, KC1, A1C1, and A1F, and the cyanide 
NaCN are among the first metal-containing molecules detected 
in space (ICernicharo & Guelinl l987; ITurner et al.1 1994). Metal- 
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bearing molecules are rarely observed in the gas phase of the 
interstellar medium due to their large refractory character that 
make them to easily form solid condensates. Their presence in 
IRC +10216 is largely caused by the high temperatures and den- 
sities prevailing close to the star, where, under thermochemical 
equilibrium, a good fraction of metals are in the gas phase form- 
ing this kind of molecules. 

The derived abundances relative to H2 in the inner layers are 
1.8 x lfr 9 for NaCl, 7 x 1(T 10 for KC1, 7 x 1(T 8 for A1C1, 1(T 8 
for A1F, and 3 x 10~ 9 for NaCN. For some of these molecules 
the abundance has been decreased in the 1-3 R, region, in or- 
der to better account for the observed line profiles. The rationale 
for this abundance modification is related to the thermochemi- 
cal equilibrium calculations, which for all these species show a 
strong abundance gradient in this region, with low values close 
to the stellar surface and higher values around the quenching re- 
gion (see Fig. 13 1. For NaCl and KC1, the abundance decrease 
in this 1-3 R* region improves significantly the agreement be- 
tween calculated and observed line profiles, indicating that these 
species are likely depleted close to the photosphere by at least a 
factor of ~ 100 for NaCl and -2-3 for KC1. For A1C1 and NaCN, 
the abundance decrease adopted in the 1-3 R* region affects only 
slightly the profiles of the observed lines, while for A1F no abun- 
dance modification has been adopted since line profiles are al- 
most not affected. The observed lines of A1C1, A1F, and NaCN 
are therefore not sensitive enough to the abundance in this in- 
nermost region, and it cannot be excluded that these molecules 
have an abundance as low as predicted by TE close to the star. 
Vibrationally excited lines should allow to better probe this re- 
gion. 

The abundances calculated under thermochemical equilib- 
rium are in reasonable agreement with our derived values, al- 
though some remarks should be made. For the metal chlorides 
NaCl and KC1, the derived abundances are higher by about one 
order of magnitude than the maximum values given by thermo- 
chemical equilibrium in the surroundings of the quenching re- 
gion (see Fig. 13 1, something that could be indicative of non- 
equilibrium processes being at work to enhance the abundances 
of NaCl and KC1. In the case of A1C1 and A1F, the derived abun- 
dances are in good agreement with the thermochemical equi- 
librium calculations, pointing towards the abundance of these 
two molecules being quenched at a radius of 2.6 R*. For NaCN, 
the analysis is not so obvious as the abundance calculated un- 
der thermochemical equilibrium shows a steep positive gradient 
that reaches the derived abundance of 3 x 10 -9 at a radius of ~6 
R», somewhat farther than the typical quenching region located 
around 3 R*. The lack of chemical kinetics data on reactions in- 
volving NaCN makes it difficult to evaluate whether such a large 
quenching radius is still possible or not. 

Metals are refractory elements, so that we expect metal- 
containing molecules to deplete from the gas phase to incorpo- 
rate into dust grains in the condensation region. Surprisingly, we 
find no evidence for such depletion except for KC1, whose abun- 
dance is found to slightly decrease from its value of 5 x 10" 10 in 
the inner layers down to 3.7 x 10" 10 at a radius of 2 x 10 15 cm. 
It therefore seems that the presence of metals in the gas phase 
is not just restricted to the warm inner circumstellar regions, 
but also to the cool outer regions where a non negligible frac- 
tion survives in the gas phase. Further evidence to this respect 
comes from the observation of other metal-bearing molecules in 
the outer layers of IRC +10216, such as MgNC (G uelin et al.l 



Table 11. Isotopic ratios in IRC +10216 



Isotopic ratio 



32 S/ 34 S 

32 s/ 33 s 

28 Si/ M Si 

28 Si/ 30 Si 
35 CI/ 37 CI 



Value 


Literature 


Solar 8 


35 ± 3.5° 


45 ±3" 


89.4 


22 ± 2.5" 


21.8 + 2.6 e 


22.1 


112 ± 12* 


121 ± 21 e 


124.9 


18 ± 2 C 


>15.4 e 


19.7 


27 ± 3.0 r 


>22.3 e 


29.9 


2.9 ± 03 d 


2.3 ± 0.24" 


3.13 




3.3 + 0.3' 





" From CS; * from SiS; c from SiO and SiS; d from NaCl, 
KC1, and A1C1; e IKahane et all (2000) an d ICernicharo eTaLl (2000) ; 
/ |Agundez etai"l (201 1); g |Asplund et aT1 (2009). 



1993), MgCN (|Ziurys et al.| 1995), A1 NC (|Ziurys et al. 
KCN (PulliametaLI 2010), and FeCN dZack et all 201: 



2002) 
), and 



from the remarkable detection of gas phase atomic metals in the 



outer envelope of IRC +10216 by |Mauron & Huggins (2010). 
As already pointed out by these latter authors, metals must be 
mostly in solid state in the outer layers of IRC +10216 since gas 
phase species do only account for a relatively small fraction of 
the total available elemental abundance, 24 % as Nal and Nail 
and 0.14 % as NaCl and NaCN in the case of sodium, 5 % as 
Kl and KII and 0.5 % as KC1 and KCN as concerns potassium, 
and 1.4 % as A1C1, A1F, and A1NC in the case of aluminium. 
These small percentages are however remarkable, given the low 
temperatures prevailing in the outer layers, and are high enough 
to allow for a metal-based gas phase chemistry to take place. In 
the inner layers the situation should not be very different, as in- 
ferred from the derived molecular abundances, with metals being 
mostly in the form of solid condensates (only beyond the con- 
densation radius) and gas phase atoms, and a small but sizable 
fraction being locked in gas phase molecules. 

5.2. Isotopic ratios 

For some of the studied molecules, a good number of transi- 
tions of rare isotopomers have been observed (see Tables [T5-[6j, 
allowing us to determine the corresponding isotopic abundance 
ratios. For each of these rare isotopomers we have adopted the 
abundance profile derived for the main isotopomer, decreased 
by the corresponding isotopic ratio, which is varied until the 
line profiles calculated with the radiative transfer model repro- 
duce the observed ones. As compared with the previous method 
used to derive isotopic ratios from line intensity ratios (e.g. 
ICernicharo et al.l 2000; Kah ane et al.l 2000), the method used 
here, based on LVG radiative transfer calculations, can better, 
although not completely, deal with optical depth effects present 
in the lines of some main isotopomers. Such effects can be com- 
pletely avoided by using optically thin lines of doubly substi- 
tuted isotopic species, but such lines are weak and are observed 
with low signal-to-noise ratios, which may introduce large un- 
certainties in the derived isotopic ratios. Here, we have therefore 
focused on single substituted isotopomers. 

The derived isotopic ratios are given in Table [TT| Under the 
reasonable assumption of a similar excitation and spatial extent 
for two isotopomers, most of the errors that usually affect the 
determination of an abundance, cancel when deriving an isotopic 
ratio (see Kaha ne et al.1 2000). The estimated uncertainties given 
in Table 1 1 therefore, neglect any error coming from the model 
and are simply based on the signal-to-noise ratio of the observed 
lines (of main and rare isotopomers), and on the sensitivity of 
the calculated line profiles to the adopted isotopic ratio. 

We find an elemental 12 C/ 13 C ratio of 35 + 3.5 from the ob- 
served lines of 12 CS and 13 CS (see Table [I]), somewhat lower 
than the value of 45 + 3 derived from intensity ratios of opti- 
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cally thin lines of SiC 2 and CS dCernicharo et al.l 2000). In the 
case of the sulfur isotope 32 S, we derive values of C 32 S/C 34 S = 
16 ± 2.5 and Si 32 S/Si 34 S = 22 ± 2.5. The different values should 
result, rather than from 32 S isotopic fractionation which is un- 
likely in warm regions, from optical depth effects that are more 
important for CS than for SiS. We thus adopt the Si 32 S/Si 34 S 
isotopic ratio as a better proxy of the elemental 32 S/ 34 S abun- 
dance ratio. A similar situation is found for the sulfur isotope 
33 S, for which we derive C 32 S/C 33 S = 83 ± 11 and Si 32 S/Si 33 S 
= 1 12 ± 12, the latter value being a better estimation of the el- 
emental 32 S/ 33 S. The isotopic ratios found for carbon and sul- 
fur are similar or slightly smaller than those previously derived 
bv lCernicharo et aLl (2000) from intensity ratios of optically thin 
lines. The fact that the 12 C/ 13 C and 32 S/ 33 S ratios derived here are 
somewhat smaller than those found by Cernich aro et al.l (2000) 
probably indicates that our results are affected to some extent 
by optical depth effects. The usefulness of the method utilized 
here is clear in cases where optically thin lines (e.g. of doubly 
substituted isotopomers) are close or below the detection limit 
of the observations. This occurs in the cases of the 28 Si/ 29 Si and 



28 Si/ 30 Si ratios, for which we find values (see Table 
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tent with the lower limits given by Cernicha ro et al. 1 (2000). 

The 35 C1/ 37 C1 isotopic ratio is found to be 2.9 + 0.3 from the 
lines of NaCl, KC1, and A1C1, which is in between the value of 
2.3 ± 0.24 found by Kah ane et al.l (2000) from intensity ratios of 
some lines of NaCl and A1C1, and the value of 3.3 + 0.3 derived 



by Agundez et al. (2011) from the modeling of the 7=1-0 
to J = 3 - 2 lines of HC1. The low 35 C1/ 37 C1 ratio found by 
IKahane et all (2000) indicates that lines of NaCl and A1C1 could 
be affected by optical depth effects, so that their calculated value 
based on line intensity ratios would be a lower limit to the true 
isotopic ratio. Finally, as concerns the 16 0/ 17 and 16 0/ 18 ra- 
tios, we have not attempted to determine them due to the low 
signal-to-noise ratios of the observed Si n O and Si ls O lines (see 
Table |2j. We nevertheless find that the observed lines of these 
two rare isotopomers are consistent with the oxygen isotopic ra- 
tios 16 0/ 17 = 840 and 16 0/ 18 = 1260 given by IKahane etaLl 
(2000). 

The isotopic ratios derived here, although systematically 
smaller, are very close to the solar values, to the exception of 
12 C/ 13 C (see Table 111. since a bias towards smaller ratios is ex- 
pected, due to the residual optical depth effects, we conclude that 
the elemental isotopic ratios 32 S/ 34 S, 32 S/ 33 S, 28 Si/ 29 Si, 28 Si/ 30 Si, 
and 35 C1/ 37 C1 in IRC +10216 are compatible with the solar val- 
ues. This result is not unexpected since the central AGB star is 
not massive enough to form much S, Si, or CI. The lower 12 C/ 13 C 
ratio found in IRC +10216, as compared with than in the Sun, 
is known to result from CNO cycle nuclear processing in the 
interior of the star (Kaha neet al.l 2000). 



6. Summary 

Astronomical observations of the carbon star envelope IRC 
+ 10216 extending over the full frequency coverage of the IRAM 
30-m telescope (80-357.5 GHz) have been analyzed to obtain 
accurate abundances for several molecules formed in the inner 
layers: CS, SiO, and SiS, the metal halides NaCl, KC1, A1C1, 
and A1F, and the metal cyanide NaCN. The observation of a large 
number of rotational transitions covering a wide range of energy 
levels, including highly excited vibrational states, allow us to 
derive abundances from the innermost layers out to the outer re- 
gions where molecules are photodissociated by interstellar ultra- 
violet photons. For some molecules, noticeably CS and SiS, we 
find that abundances are depleted from the inner to the outer lay- 



ers, implying that they contribute to the formation of dust. The 
amount of sulfur and silicon locked up in gas phase molecules 
in the inner layers, 27 % and 5.6 % respectively, indicates that a 
major fraction of these elements is in the form of solid conden- 
sates, most likely as MgS and SiC, in IRC +10216's envelope. 
Metal-bearing molecules lock a relatively small fraction (a few 
percent) of metals. Nevertheless, we find that NaCl, KC1, A1C1, 
A1F, and NaCN, which all contain refractory elements, seem to 
maintain a surprisingly high abundance in the cold outer layers, 
where they can participate, together with metallic atoms, in a 
rich gas phase chemistry. The lines of rare isotopomers of CS, 
SiO, SiS, and some metal compounds have been analyzed, and 
we derive abundance ratios for rare isotopes of carbon, sulfur, 
silicon, and chlorine. The values obtained are in good agreement 
with previously reported values. They are consistent with the so- 
lar values, except for the 12 C/ 13 C ratio, which is more than twice 
smaller in IRC +10216 than in the Sun, in agreement with pre- 
vious findings. 
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Fig. 11. Rotational lines of NaCN in IRC +10216 as observed with the IRAM 30-m telescope (black histograms) and as calculated 
with the radiative transfer model (blue lines). Lines are ordered by increasing frequency from bottom to top (and from left to right). 
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Table 1. CS line parameters in IRC +10216 



Transition 


Vrest 


Vobs 


JT*dv 


^exp 


Notes 




(MHz) 


(MHz) 


Kkms- 1 


km s 








CS 








v=0 J=2-l 


97980.952 


97980.8(1) 


117(3) 


13.8(2) 




v=0 J=3-2 


146969.025 


146968.6(2) 


267(6) 


13.8(3) 




v =0 J=4-3 


195954.213 


195953.7(4) 


286(8) 


13.3(3) 




v=0 J=5-4 


244935.554 


244935.1(3) 


308(8) 


12.5(4) 




v=0 J=6-5 


293912.089 


293911.7(3) 


423(10) 


12.6(4) 




v=0 J=7-6 


342882.854 


342881.8(3) 


259(8) 


12.3(4) 




v=l J=2-l 


97270.997 


97271.5(5) 


0.11(2) 


6.1(6) 


$! 


v=l J=3-2 


145904.089 


145903.8(3) 


1.05(5) 


5.4(4) 




v=\ J=5-4 


243160.648 


243161.1(5) 


1.38(7) 


6.6(3) 


S 


v=l J=6-5 


291782.190 


291782.2(5) 


3.0(1) 


11.8(8) 




v =l J=7-6 


340397.957 


340398.3(5) 


1.74(8) 


10.4(6) 




v=2 J=6-5 


289651.545 


289651.9(5) 


0.60(5) 


5.8(6) 




v=2 J=7-6 


337912.189 


337912.9(6) 


0.34(3) 


4.8(6) 




v=3 J=5-4 


239608.923 


239610.1(5) 


0.07(1) 


3.3(2) 




v=3 J=6-5 


287520.096 


287521.1(8) 


0.28(2) 


7.1(5) 




v=3 J=7-6 


335425.481 


335425.8(8) 


0.14(2) 


3.4(5) 




v=0 J=2-l 


92494.273 


"CS 
92494.2(2) 


4.1(2) 


14.2(8) 


S 


v=0 J=3-2 


138739.267 


138739.1(2) 


10.7(5) 


14.1(4) 


s 


v=0 J=5-4 


231220.689 


231220.6(2) 


16.8(8) 


14.2(3) 




v=0 J=6-5 


277455.403 


277456.2(6) 


22.6(10) 


14.5(6) 




v=0 J=7-6 


323684.976 


323686.3(10) 


12.5(8) 


14.4(8) 




v=l J=5-4 


229592.986 


229593.0(6) 


0.05(1) 


2.5(2) 




v=0 J=2-l 


96412.956 


C J4 S 
96412.9(2) 


10.1(4) 


14.3(3) 




v=0 J=3-2 


144617.107 


144617.0(1) 


21.5(5) 


13.8(2) 




v=0 J=5-4 


241016.098 


241016.1(1) 


39.4(8) 


14.4(2) 




v=0 J=6-5 


289209.077 


289208.9(5) 


40(1) 


13.7(4) 




v=0 J=7-6 


337396.471 


337396.1(5) 


24.2(8) 


14.0(4) 




v=l J=3-2 


143577.671 


143576.9(10) 


0.04(1) 


6.9(10) 


1 


v=l J=6-5 


287130.180 


287131.1(5) 


0.18(2) 


4.9(7) 





v=0 J=2-l 
v=0 J=3-2 
v=0 J=5-4 
v=0 J=6-5 
v=0 J=7-6 



97172.063 
145755.730 
242913.608 
291485.927 
340052.572 



C jj S 
97171.6(5 
145755.6(4 
242913.7(2 
291485.9(6 
340052.8(8 



1.71(8) 
4.9(2) 
8.3(1) 
10.0(3) 
6.1(2) 



15.1(7) 
14.4(7) 
14.4(4) 
14.2(6) 
13.6(8) 



Notes - Numbers in parentheses are lcr errors in units of the last digits. 
Observed frequencies v obs were derived adopting a systemic velocity v LSR for 
IRC +10216 of -26.5 km s~» dCernicharo etaLl 2000). 
|| Line width parameter v cxp has been fixed. 

$ Blended with another line. An individual fit for each line is possible. 
£ Two transitions are very close in frequency. A single line is fitted. 
! Marginal detection. 
* Clear or probable maser line. 

oo Complex line profile. Line fit is not particularly good. 
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Table 2. SiO line parameters in IRC +10216 



Transition 


y rcst 




jTJdv 


v cxp Notes 




(MHz) 


(MHz) 


K km s" 1 


km s -1 






SiO 






v=0 J=2-l 


86846.971 


86846.8(2) 


52(2) 


13.9(2) 


v=0 J=3-2 


130268.665 


130268.5(2) 


106(4) 


13.6(2) 


v=0 J=4-3 


173688.210 


173687.6(3) 


153(6) 


13.6(3) 


v=0 J=5-4 


217104.889 


217104.6(2) 


145(7) 


13.1(5) 


v=0 J=6-5 


260517.985 


260517.8(2) 


161(8) 


13.7(4) 


v=0 J=7-6 


303926.783 


303926.6(2) 


173(8) 


13.0(4) 


v=0 J=8-7 


347330.565 


347329.9(7) 


153(10) 


13.0(7) 


v=l J=2-l 


86243.430 


86243.1(10) 


0.04(2) 


10.8(10) $ ! 


v=\ J=5-4 


215596.029 


215596.0(5) 


0.54(5) 


13.4(8) 


v=l J=7-6 


301814.369 


301815.9(10) 


0.93(6) 


12.6(8) 






2lJ SiO 






v=0 J=2-l 


85759.202 


85759.0(2) 


4.2(2) 


15.3(4) 


v=0 J=4-3 


171512.814 


171513.0(2) 


11.9(5) 


14.1(3) 


v=0 J=5-4 


214385.778 


214385.9(3) 


13.2(5) 


13.9(4) 


v =0 J=6-5 


257255.249 


257255.7(5) 


11.1(5) 


14.0(5) 


v=0 J=7-6 


300120.528 


300121.7(3) 


16.4(6) 


13.8(4) 


v=0 J=8-7 


342980.917 


342980.9(5) 


11.2(5) 


13.4(5) 






3U SiO 






v=0 J=2-l 


84746.193 


84746.0(2) 


2.9(1) 


14.6(3) 


v=0 J=4-3 


169486.930 


169486.8(2) 


7.3(3) 


14.5(3) 


v=0 J=5-4 


211853.546 


211853.4(3) 


15.2(6) 


13.5(5) 


v=0 J=6-5 


254216.751 


254216.7(6) 


9.5(5) 


14.6(6) 


v=0 J=7-6 


296575.863 


296575.9(4) 


10.4(2) 


14.0(5) 


v=0 J=8-7 


338930.199 


338930.0(4) 


6.9(1) 


14.0(5) 






Si 1!i O 






v=0 J=2-l 


80704.907 


80705.6(15) 


0.06(3) 


14.5(20) ! 


v=0 J=5-4 


201751.443 


201750.8(10) 


0.30(8) 


12.6(15) ! 






Si lv O 






v=0 J=4-3 


167171.974 


167172.1(6) 


0.49(3) 


14.3(6) 


v=0 J=5-4 


208959.990 


208960.4(8) 


0.30(5) 


14.2(6) 


v=0 J=7-6 


292525.404 


292526.0(10) 


0.54(3) 


14.1(8) 


v=0 J=8-7 


334301.474 


334300.7(10) 


0.33(3) 


16.4(20) 
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Table 3. SiS line parameters in IRC +10216 



Transition 


Vrest 


Vobs 


1 T* dv 


v exp 


Notes 




^ivi nz j 


\1\LlTZ,} 
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km s 








515 








lf _n I — S A 
V— U J —3-4 


Qfi771 SAA 
yU / / 1 .304 
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y\J III .3(1 ) 
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1 A 

14.U(,3 ) 




ft r_ A c 

V— U J— OO 


1 0SQ7A 301 


1 OSQ7A r \(\ \ 
L\JoyZ £ +.j\l ) 


oypj 


1 A 1 (A\ 
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lf — ft T— 8 7 
V— U J — o- / 
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1 S1 S7S C\(l\ 
lol3Z3.U^Z^ 


1 7S/'S , i 
1 / OO ) 


1 7 Qr7A 




ai— n T— 1 1 1 n 
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7 1 78 1 7 AS3 

Zl / 1 / .033 


7 1 70 1 7 ^;('7^ 
Zl /ol /.3^Z^ 


1 78^^ 
1 / 0^0 ) 


13. "4^0 ^ 




i,_ ft T— i a n 
V— U J — 13-1Z 


73SQA1 3A3 
Z33V01 .303 


7^SQAft Ql'^^ 
ZjjyOV.y{j ) 


7 1 S/'AA 
Z13^0 J 


1 7 Q/'SA 

iz.y^o ^ 




ai— n T— 1 /i n 

V— U J — 14-1 J 


7S/1 1 fl3 7 1 1 
Z341U3.Z1 1 


7S/i 1 m c\( 1 n^i 
Z341U3.U(,1U) 


oni / 1 7\ 
3U1(,1Z) 


11 1 / 1 c\ 

13.1(13) 


CO 4* 


a,— n i— 1 s i a 

V— U J — 1 J- 14 


7777A3 flS7 
Z / ZZ43.U3Z 


7777A7 f\( 1 O^i 
z / zz^z.OyiV) 


Q7 1 / 1 a\ 
3 / 1(_14J 


17 S^l SA 
1Z.3( 13 ) 


CO A 


a,— n T— 1 A is 
V— U J — 10-1 3 


700380 7AA 
ZyUJoV. /44 


7on^sn ka\ 

Zy\) 3oU. /{,4_J 


Q^^/ 1 AA 


1 2 S("AA 
13.3(0) 




ai— n T— n i a 

V— U J — 1 /-10 


•2fl8S1A 1 AA 
3U03 10. 144 


QfiSS 1 S ^f'A 1 ! 
3Uo3 13.31^0 ) 


1 SS/'SA 
103^0 J 


1 7 Q('7A 




ai— n T— 1 s 1 1 
v—y) J — lo-l / 


37AA/1Q 1 HS 
3Zu04V. 1U0 


3Z0045. 1 yo) 


707/1 Q\ 

zjz(lj) 


1 1 f\(Q\ 
13.U(o) 




a,— n T— 1 Q is 

v— u j — ly-io 


3AA77Q AQ7 
344 / fy.'-tyZ 


^AA77S Q{&\ 
344 / / o.7^0 ) 


1 QA/^ 1 1 A 
1V4^1 1 ) 


13.Z(0) 




Ai—1 T— S /I 
V— 1 J— j-4 


yu3zy.oyu 


yvjjzy.yyZ) 


KJ.zyyZ) 


1 1 7("AA 
ll./(0) 




V— i J — no 


1 OS3QA 701 

iuo3y4.zy i 


1 ns^QA 7^1 ^ 
!Uo3y4. /^l^ 


H AS/^AA 
U.00(,4; 


Q 7/"AA 

y. 1(0) 




Al— 1 T— G 7 

V— 1 J— o- / 


1 A A S70 3A7 
1443ZU.30 / 


1 A A S7fl S^/I^i 
1443 ZU. 3(4 ^ 


U.o3(,3 ) 


A S/AA 
O.3(0J 




,._ 1 T_Q O 

V— 1 J— y-o 


1A7SS1 7SS 
10Z30 1. /33 


1 A7SSO A^l f\\ 
10Z3oU.4^1U^ 


1 7Ai'inA 
1 .Z4^1UJ 


Q Sl( 1 HA 




,._ i t— 1 1 1 ft 

V— 1 J — 1 1-1U 


1 08700 S 1 7 

iyo /uu.3 1 / 


1 QS7fin 1(A\ 
lyo /UU.Z^4^ 


1 A7^SA 

1 .ozp ) 


1 7 (VSA 
1Z.U(3 ) 




V— 1 J — 1Z-1 1 


7 1 A7S7 A03 
Z10 /3 / .0U3 


7 1 A7S7 KA\ 
Z10/3 / . /^4^ 


7 HS^AA 
Z.UO(,0 ) 


1 7 Qf'AA 
1Z.V(0) 




Al— 1 T— 1^1 17 

V— l J — 13-1Z 


73/1 8 1 7 OAS 
Z3401Z.V0O 


Z34ol3.U(3) 


7 AA^CA 
Z.00(,0 ; 


1 fl 7("SA 
1U. /(3) 




Al— 1 T— 1A 1^ 

V— 1 J — 14-1 J 


7S78AA AAS 
Z3Z000.400 


7S7SAS 1 (\ (Y\ 
Z3Zo0o.H,lU^ 


7 Q^i' 1 OA 
Z.y3^1UJ 


17 ^^1 HA 
1Z.3( 1U) 




il— 1 I— 1 ^ 1 A 
V— 1 J — 13-14 


77fiO 1 7 OA 1 
Z /UV 1 /.V01 


77flO 1 A Q^'7^ 
Z IVylO.yy 1 ) 


3.3(1) 


1 ft /I /"SA 
1U.4(0 ) 




Al— 1 T— 1 A 1 S 

V— 1 J — 10-1 J 


7SSQA7 ^f!7 
ZooVO / .3UZ 


7SSQA7 Sl'A^ 


A Sf7A 
0.3(Z) 


1 7 Tf'AA 
1Z.Z(0) 




ai— 1 T— 1 Q 1 S 

v— i j — iy-io 


•I -i q i nn qsa 

343 1UU.704 


343 IUU.4^0 ) 


^ 7^ 1 A 
5.1(1) 


1 1 ^^7A 
11. 5(1) 




u—7 I — S A 
V— Z J —3-4 


SQSSS 1 7S 
07000. 1/3 


SQSSS ACX\ 


n no/^ 1 a 


Q AfSA 




At— 7 I— A ^ 

V— Z .1—00 


1 H7SAA 7^ 1 
1U / 004. Z3 1 


1 H7SAA 1 (1\ 
1U / 004. 1 \Z) 


KJ.iyyZ) 


1 ft Afl\ 
1U.4( / ) 




it— 7 1—8 7 
V— Z J — 5- / 


1 A^S 1 ^ A1 S 
1430 13.013 


1 A^S 1 ^ (VS^i 
1430 13.U^3 ^ 


n 1 A^7A 


S Q/'SA 




ii—7 i— 1 1 1 n 

V— Z J — 1 1-1U 


1 Q777S 7 1 S 

iy / /zs. / lo 


1 Q777S K£\ 
ly 1 1 zo. 1(0) 


U.3 /(,4^ 


1 ft A (Q\ 




11—7 T— 1 7 11 
V— Z J — 1Z-1 1 


7 1 SAQ7 AS7 
Z130V / .43Z 


7 1 SAQ7 Af^ 
Z1307 / .4^3 ) 


n SA/^AA 
U. 34(^4 y 


A 7fAA 
O. 1(0) 




u—7 i— i Q n 

V— Z J — 1 3-1Z 


7^^AAA AAA 
Z33004.404 


7^^AA^ f\( 1 C\\ 


n A/VAA 


A 1 1 HA 
0. 1(1U) 


<C 
3> 


ii—7 I— 1/1 1 
V— Z J — 14-1 J 


7S 1 A7Q AOQ 

Z3 lozy.ouy 


7S 1 A7S Ai'S^ 
Z3 lOZo.O^o ) 


1 A7^SA 
1 .hZ(o ) 


Q AfSA 

y.o(t$) 




ii—7 I— 1 S 1 A 
V— Z J — 13-14 


7AQSQ7 7AA 
Z0V3yZ. /44 


7AQSQ^ A(9\ 
Z0V3y3.4^o ) 


1 SQ^7A 


7 AfSA 




ii—7 T— 1 A 1 S 
V— Z J — 10-1 J 


7S7SS^ 77A 
Zo /333. /ZO 


7S7ssi Anru 

ZO /331.0(,1U^ 


U. 44(^4 7 


S f\( 1 HA 
O.O(IU) 




ii—7 T— 1 7 1 A 
V— Z J — 1 /-10 


•2nSS 1 7 A 1 fl 
3U33 1Z.41U 


anSS 1 7 
3U33 1Z. I\l ) 


1 lA^^A 
1.14(J) 


A Sf7A 




it—7 T— 1 Q IS 

v— z j — iy-io 


^A1 A77 ^1 A 
3414ZZ.3 14 


3414ZZ.4^ / ) 


n s 1 /^a 

U.3 YyJ ) 


S ftfSA 




ai— ^ 1— S A 
V—J J —3-4 


SQAAA AHA 
O7440.4U4 


SQAAA ^(A\ 
07440.3^4^ 


U.UZ3^ / ) 


^ S/^1 HA 
3.3(1U) 




ii-^l I— A S 
V — 3 J— 0-3 


1 H7^^A 1 flA 
1U /334. 1U4 


LXJI JJJ.\J(L\J) 


\J.Z\jyj ) 


s rtf 1 n^i 




n—^l T— 8 7 


i /i q i r\c. nn a 
143 1U0. / /4 


1 A 1 1 H7 

143 1U 1 .j(o ) 


n 1 1 C7A 

U. 1 lyz) 


A 1 /SA 

4.1(8) 




ii—^ T— Q S 
V—J .1 — V-o 


1 AOQQ1 ASA 
IOU77 1.430 


1AOQQ1 ^i'S^ 

iouyyi.3^o^ 


n 1 S^7A 

u. ijyZ) 


A 1 fSA 
4.1(5) 




li—^i T— 1 n q 
v—j j — lu-y 


1 7SS7A SAS 
1 / oo /4.o43 


1 7SS7S Ai'S^ 
1 / 00 /3.4(^o ) 


n 1 7^7A 

u. izyz ) 


^ 7f7A 
3. /(/) 




n— ^ T— 1 1 1 ft 
V—J J — 1 1-1U 


1 QA7SA 7QA 
1^0 / 30. fyO 


1 QA7S7 C\< 1 

iyo /3 / .u(,iu^ 


n ns^7A 


7 SfAA 
Z.3(0) 




il— 1 T— 1 7 11 
V—J J — 1Z-1 1 


7 1 A AQ7 1 A7 
Z1403 / . 10 / 


7 1 /I ATA /I (Q\ 
Z1403u.4(o) 


n 1 S/7A 

U.lo(,zj 


7 7/AA 




a T i a 1 1 

V—J J — 14-13 


7Sft^Q7 SQ1 
ZjyjJyZ.jy 1 


7so^qa ^^'7^^ 

Z3U3V4.U(,ZU^ 


U.3 /(,0; 


S ftfSA 
O.U(0 ) 




ii— T T— 1 S 1 /I 
V— J J — 13-14 


7AS7A7 ^1S7 
ZOozO /.33 / 


7AS7A7 S/ 1 

Z0oz0/.3(,1U) 


H 7S/^/1A 


A A/1flA 
4.0(1U) 




ii—^ T— 1 A i 

V— 3 J — 10-13 


7SA1 QA7 
ZS0137.70 / 


7SA1 AO AfH\ 
Zo014U.4^o ) 


n ^S/^AA 
U. 30(^4; 


S ftfSA 
3.U(0 ) 




ii— 'i T— 1 7 1 A 
V— D J — 1 /-10 


an/in 1 n 77S 
3U4U1U.Z/J5 


3U4U1U. i(lz) 


U. /3(,o; 


S 1< 1 7A 
3.Z(1Z) 


5> 


n— a T— 1 S 17 
V— 3 J — lo-l / 


■37 1 070 1 Af\ 
3Z1 / 0. 140 


Q7 1 070 7/ 1 C\ 
3Z10 lo. 1(1 J) 


n AA/^SA 
U.40O ) 


A 0/ 1 SA 
4.y(13 ) 


3> 


A 1—8 7 
V— 4 J— o- / 


1 A7^QQ S 1 Q 

i4Z3yy.oiy 


1 A7AHO 7^1 O^i 
1 4Z4UU . Z^ 1 


n ns/^ 1 a 

U.U3(,1 ) 


^ 7fSA 
3.2(8) 




ai— 4 T— 1 7 11 
V — 1 J — 1 Z- 1 1 


7 1 ^S7A 710 
Z 133 / U. / 1U 


7 1 ^S7A Qf'S 1 ! 
Z 133 / U."^3 ) 


n 7 1 m.a 

u.z l yn ) 


^ SfSA 
3.3(3 ) 




n— A T— 1 ^17 
V— 4 J — 13-1Z 


7Q 1 ^AA Q7A 
Z3 1300.7 /O 


7a 1 •2A7 7*"^ 
Z3 130 / .Zl^3 ) 


n 7A/^A 
U.Z0(,3 ; 


7 7f^A 

z. /(3) 




A T— 1 A 13 
v — *+ J — 1H— 1 3 


7AQ1 SS ^77 
Z^+y lJD.j I Z 


7AQ1 SS Af'Ai 
Z"4!7l33.0^4^ 


n 17/^zlA 


^ SfAA 
3.3(4) 




ai— 4 T— 1 S 1 A 
V— 4 J — 13-14 


7AAQA1 7SA 
Z00V41. /34 


7AAQA7 Af7^ 
Z00V4Z.4^ / ^ 


U.Z4(^4; 


A ftf7A 
4.U(/) 




,,—A T— 1 A IS 
V— 4 J — 1D-13 


7SJ.77S Q7S 
Zo4 / Zj.y I 


7SJ.77S 1 
Zo4 / Z3. 1 l^o^ 


n 1 An a 
u. lot, 17 


^ 7fQA 
3.Z(V) 




ai A T— 1 7 1 A 

V— 4 J — 1 /-10 


3UZ3U / .7Ul 


j\JZj\Jy.J(LD ) 


n 7S^AA 
U.Z0(,4; 


s n/1 sa 

3.U(13 ) 




v=4 J=18-17 


3202S7 377 


320286 4f20i 


0.09(2) 


3.3(10) 


I 


v =4 J=19-18 


338064.264 


338066.3(20) 


0.15(3) 


5.0(15) 




v=4 J=20-19 


355838.417 


355839.9(15) 


0.22(3) 


7.3(15) 




v=5 J=14-13 


247917.909 


247917.2(15) 


0.16(2) 


6.4(20) 




v=5 J=16-15 


283311.709 


283313.1(20) 


0.16(2) 


2.9(10) 




v=5 J=17-16 


301005.224 


301005.5(20) 


0.17(4) 


5.5(20) 


s 


v=5 J=18-17 


318696.291 


318697.2(15) 


0.13(3) 


4.7(20) 


$ 


v=5 J=19-18 


336384.765 


336385.1(8) 


0.12(2) 


4.5(10) 
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Table 3. Continued. 



Transition 


v rest 


y obs 


]'T A dv 


v^p Notes 








(MHz) 


(MHz) 


Kkms" 1 


km s _1 










2<) SiS 






V 


-0 T- 
— U J — 




89103.780 


89103.8(1) 


3 9(1) 


14.8(5) 


V 


-0 T- 
-U J — 


00 


106923 018 


106923 0(1) 


1 8(1) 


15.1(5) 


V 


-0 T- 
— u J — 


8 7 
o- / 


142558.872 


142558.7(2) 


14.3(3) 


14.6(3) 


V 


— U J — 


s 
7-0 


160375 212 


160375 1(1) 


19 0(4) 


14.0(3) 


V 


-0 T— 
— u J — 


11 10 

1 1- 1U 


196004.027 


196004.1(2) 


16.6(4) 


14 1(3) 


V 


-0 T— 

— U J — 


19 11 
1 L- 1 1 


213816.227 


213816.1(2) 


20.1(4) 


14.0(3) 


V 


-0 I- 

— U J — 


1 ^ 19 
1 J- LZ, 


231626.770 


231626.7(2) 


24.0(7) 


13.8(3) 


V 


-0 T— 
— u J — 


u ia 
i L +- 1 j 


249435 519 


249435.5(3) 


20.6(5) 


13 4(4) 


V 


-0 T— 
— u J — 


1 J-lf 


267242.336 


267242.0(4) 


23.2(6) 


14.6(5) 


V 


-0 T— 

— U J — 


1 A 1 S 
10-1 J 


285047 083 

— 1 ' ~ ' V./ 1 / . V./ "..1 ^ / 


285046 8(5) 


24.9(10) 


13.7(6) 


V 


-0 T- 

— U J — 


17 16 
1 /-io 


302849.622 


302850.3(7) 


22.8(8) 


13 9(7) 


V 


-0 T- 
— U J — 


18 17 
1 o-l / 


320649.814 


320649.1(5) 


20.0(7) 


13.9(6) 


V 


-0 T- 
— U J — 


1 Q 1 8 
17-10 


338447.522 


338447.5(5) 


11 7(5) 


14 2(5) 


V 


- 1 T— 
— 1 J — 


Q 8 
y-o 


159602 038 


159600 6(15) 


0.12(3) 


8.1(10) ! 


V 


- 1 T— 
— 1 J — 


M n 

1H— 1 J 


248232.767 


248232.9(5) 


0.27(4) 


6 1(6) 




= 1 J= 


16-15 


283672.487 


283673.7(8) 


0.38(2) 


5.6(6) 


V 


= 1 J= 


17-16 


301389.100 


301390.1(6) 


0.46(3) 


6.4(7) 




= 1 J= 


18-17 


319103.365 


319105.5(20) 


0.50(6) 


7.5(10) 


y 


= 1 J= 


19-18 


336815.143 


336816.8(15) 


0.40(4) 


10.1(10) 




=2 J= 


12-11 


211754.292 


211755.1(8) 


0.07(2) 


4.4(7) ! 




=2 J= 


13-12 


229392.993 


229392.4(8) 


0.05(1) 


1.4(3) ! 


V 


-9 T- 

— Z, J — 


1 6-1 S 


282297.754 


282298.9(15) 


0.08(1) 


2.6(8) $ 




-2 J= 


17-16 


299928.433 


299928.8(10) 


0.07(2) 


2.2(5) 




=2 J= 


18-17 


317556.761 


317556.0(10) 


0.17(3) 


7.6(10) 


V 


-0 T— 
— U J — 


S zl 


87550.623 


3U SiS 
87550.5(3) 


2 5(1) 


14.3(5) 


V 


-0 T— 
— U J — 


00 


105059.283 


105059.2(2) 


4 9(1) 


14 2(3) 


V 


-0 T— 
— U J — 


8 7 
o- / 


140074 069 


140074 1(4) 


7 7(2) 


14 7(5) 


V 


-0 T- 
— U J — 


Q 8 
y-o 


157579.931 


157579.8(3) 


12.2(3) 


14 1(3) 


V 


-0 T- 
— U J — 


10 Q 


175084 593 


175084 1(3) 


7 4(2) 


14.3(3) 


V 


-0 T- 
— U J — 


19 11 

1 L- 1 1 


210089.787 


210089.6(3) 


14 4(4) 


14 2(3) 


V 


-0 T— 
— U J — 


1 1 19 


227590.052 


227589.9(3) 


15.7(5) 


13.6(5) 


V 


-0 T— 
— U J — 


1 ZL 1 ^ 
1H— 1 J 


245088.585 


245088.7(3) 


14.3(5) 


13.9(3) 


V 


-0 T— 

— U J — 


1 S 1 A 
1 J-lf 


262585 253 


262585.0(3) 


13.5(6) 


14.3(3) 


V 


-0 T- 

— U J — 


1 f\ 1 S 

1 U 1 J 


280079.923 


280080.8(6) 


16.2(8) 


14.0(5) 


V 


-0 T— 
— U J — 


17 1 fx 
1 /-IO 


297572.461 


297572.1(5) 


13.2(6) 


14 0(4) 


V 


-0 T— 

— U J — 


18 17 
l o- 1 / 


315062.733 


315062.1(5) 


11.8(5) 


13.5(5) 


V 


-0 T— 
— U J — 


1 Q 1 8 
ly- Lo 


332550 608 


332549 9(7) 


10.6(6) 


14 2(5) 


V 


-0 T- 

— U J — 


90 1Q 

ZU 1 7 


350035.950 


350035.7(4) 


8 1(3) 


13.8(5) 


V 


-1 T- 
— 1 J — 


8 7 
o- / 


139404.721 


139404.4(8) 


0.05(1) 


5 1(8) 


V 


-1 T— 

— 1 J — 


1^19 


226502.331 


226502.3(3) 


0.20(2) 


5 4(6) 


y 


= 1 J= 


14-13 


243917 185 


243916 7(5) 


0.24(2) 


6.6(5) 


V 


= 1 J= 


16-15 


278741.159 


278741.3(10) 


0.42(5) 


8.8(10) $ 


V 


= 1 J= 


17-16 


296150.011 


296149.2(8) 


0.20(2) 


3.2(7) 


V 


= 1 J= 


19-18 


330960.781 


330960.3(6) 


0.31(3) 


3.8(8) 


V 


=2 J= 


12-11 


208081.602 


208083.1(10) 


0.06(1) 


6.4(8) 


V 


=2 J= 


14-13 


242745.671 


242745.2(5) 


0.09(1) 


3.2(5) 


V 


=2 J= 


17-16 


294727.422 


294727.9(7) 


0.09(1) 


3.9(6) 


V 


=3 J= 


13-12 


224326.538 


224327.2(5) 


0.03(1) 


3.2(4) ! 


V 


=3 J= 


14-13 


241574.006 


241573.7(6) 


0.05(1) 


4.4(5) ! 



Agundez et al: Molecular abundances in the inner layers of IRC +10216, Online Material p 5 



Table 3. Continued. 



Transition 


v rest 


y obs 


]'T A d\ 




Notes 




(MHz) 


(MHz) 


K km s _1 


km s" 1 








Si J4 S 








v=0 J=5-4 


88285.834 


88285.6(3) 


2.71(8) 


14 9(4) 


$ 


v=0 J=6-5 


105941 510 


105941 5(2) 


4 7(1) 

-t. /V, 1 / 


14 2(4) 




v=0 J=8-7 


141250.286 


141250 1(3) 


9 7(2) 


13 9(4) 




v=0 J=9-8 


158903 116 


158903 0(3) 

x ^j kj y \j *j .\j y ~j j 


16 4(4) 


14 0(4) 




v=0 J=10-9 


176554.726 


176554 5(3) 


14.5(3) 


13.9(5) 




v=0 J= 12-1 1 


211853.744 


211853 5(3) 


15.4(3) 


14 1(4) 




v=0 J=13-12 


229500.881 


229501.0(4) 


15.5(5) 


13 9(7) 


$ 


v=0 J=14-13 


247146.257 


247146.1(3) 


22.2(8) 


14 2(4) 




v=0 J=15-14 


264789.734 


264789.9(3) 


19.5(5) 


13.8(4) 




v=0 J= 16-15 


282431 179 


282432.1(7) 


22.2(6) 


13.3(5) 




v=0 J=17-16 


300070.455 


300071.1(7) 


20.1(6) 


13.3(5) 




v=0 J=18-17 


317707.428 


317707.1(6) 


18.5(6) 


14.0(5) 




v=0 J=19-18 


335341.960 


335341.8(6) 


12 7(4) 


14 1(5) 




v=0 J=20-19 


352973.918 


352975.1(8) 


13.3(6) 


13 1(7) 

x .J . X ^ / ) 




v=l J=9-8 


158140.565 


158140.3(7) 


0.11(2) 


5.6(6) 




v=l J=12-ll 


210836 992 


210836.6(8) 


0.11(2) 


6.8(6) 




v=l J=14-13 


245960.030 


245960.3(4) 


0.13(2) 


5 9(6) 




v=l J=15-14 


263518.767 


263518.1(6) 


0.51(4) 


10.8(8) 




v=l J=16-15 


281075.470 


281076.4(8) 


0.34(3) 


7.2(10) 




v=l J=17-16 


298630.002 


298631.0(15) 


0.49(6) 


6.8(15) 


$ 


v=l J=18-17 


316182.227 


316184.0(15) 


0.55(6) 


8.6(15) 


S 


v=l J=19-18 


333732.011 


333732.2(8) 


0.20(3) 


5.7(8) 




v=l J=20-19 


351279.217 


351278.8(6) 


0.35(3) 


9.6(6) 




v=2 J=18-17 


314656.884 


314657.6(10) 


0.16(2) 


6.3(8) 




v=2 J=19-18 


332121.910 


332121.2(10) 


0.07(2) 


3.7(10) 


! 


v=3 J=14-13 


243587.208 


243586.9(10) 


0.05(1) 


3.0(7) 


! 


v=0 J=5-4 


89489 232 


Si 33 S 
89489 2(3) 


0.75(5) 


14 1(5) 


S 


v=0 J=8-7 


143175.531 


143175.5(5) 


2 4(1) 


14.5 


II $ 
II ■+* 


v=0 J=9-8 


161068 923 


161069 0(3) 


2.48(7) 


14.0(3) 




v=0 J=10-9 


178961.062 


178961.0(5) 


3.36(9) 


14 1(4) 


$ 


v=0 J=l 1-10 


196851 808 

J. y u u j x . vj v/ u 


196851 9f5^ 

x y vj vj *j x . y \ *j j 


2.05(8) 


14.6(5) 




v=0 J=12-ll 


214741.024 


214741.0(4) 


2.27(7) 


14 5(4) 




v=0 J= 13-12 


232628.569 


232628.4(3) 


3 5(1) 


14 U3) 




v=0 J=14-13 


250514.304 


250514.4(4) 


1.71(7) 


14.1(4) 




v=0 J=15-14 


268398.090 


268398.7(6) 


5.0(4) 


14.5(4) 




v=0 J=16-15 


286279.788 


276279.3(8) 


3.1(2) 


13.8(8) 


$ 


v=0 J=17-16 


304159.258 


304159.1(4) 


4.8(2) 


13.2(6) 




v=0 J=18-17 


322036.362 


322036.6(8) 


3.8(2) 


13.9(7) 




v=0 J=19-18 


339910.961 


339910.8(7) 


2.52(8) 


13.7(6) 




v=l J=17-16 


302689.236 


302690.2(10) 


0.18(3) 


4.6(8) 


$ 


v=l J=18-17 


320479.853 


320479.5(10) 


0.05(1) 


1.9(8) 


! 
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Table 4. NaCl line parameters in IRC +10216 



Transition 


v rcst 


y obs 




Vexp 


Notes 




(MHz) 


(MHz) 


Kkms-' 


km s" 1 








Na 35 Cl 








v=0 J=7-6 


91169.875 


91170.0(3) 


1.12(5) 


15.3(5) 




v=0 J=8-7 


104189.660 


104189.6(4) 


1.12(7) 


13.7(7) 




v=0J=10-9 


130223.627 


130223.5(3) 


1.78(8) 


13.8(6) 




v=0 1=11-10 


143237.360 


143237.4(3) 


1.95(8) 


13.1(6) 




v=0 1=12-11 


156248.627 


156248.6(4) 


2.26(10) 


13.1(7) 




v=0 J=13-12 


169257.205 


169260.1(20) 


1.24(20) 


10.8(15) 


j 


v=0 1=16-15 


208264.556 


208264.6(3) 


2.18(10) 


14.3(8) 


$ 


v=0 J=17-16 


221260.132 


221256.0(3) 


2.44(8) 


14.2(4) 




v=0 1=18-17 


234251.898 


234251.2(4) 


2.14(8) 


14.8(4) 




v=0 J=19-18 


247239.628 


247239.4(4) 


2.63(6) 


13.3(5) 




v=0 1=20-19 


260223.098 


260223.8(8) 


4.5(2) 


14.5(7) 




v=0 1=21-20 


273202.085 


273204.1(15) 


2.44(10) 


13.2(8) 




v=0 1=22-21 


286176.364 


286176.2(8) 


2.03(8) 


14.3(7) 




v=0 1=23-22 


299145.710 


299145.8(6) 


1.90(6) 


15.9(10) 




v=0 1=24-23 


312109.900 


312109.8(6) 


1.39(5) 


12.2(6) 




v=0 J=26-25 


338021.912 


338021.6(7) 


1.13(4) 


13.2(7) 




v=0 1=27-26 


350969.286 


350968.9(7) 


1.14(4) 


13.1(6) 








Na J 'Cl 








v=0 J=7-6 


89220.115 


89219.9(5) 


0.39(3) 


15.7(10) 




v=0 J=8-7 


101961.552 


101961.2(4) 


0.41(2) 


13.6(6) 




v=0 J=9-8 


114701.273 


114701.1(8) 


0.43(5) 


13.0(15) 


$ 


v=0 1=16-15 


203813.203 


203812.0(15) 


0.55(7) 


14.0(10) 




v=0 1=20-19 


254663.465 


254662.8(15) 


0.74(10) 


11.9(20) 


$ 


v=0 1=21-20 


267365.833 


267365.9(8) 


0.83(7) 


13.3(8) 




v=0 1=23-22 


292756.828 


292756.1(7) 


0.89(6) 


12.5(7) 




v=0 1=24-23 


305445.026 


305444.3(8) 


1.16(7) 


12.4(6) 




v=0 1=25-24 


318128.071 


318127.4(10) 


0.43(5) 


15.0(10) 




v=0 1=26-25 


330805.749 


330805.5(10) 


0.57(5) 


10.9(10) 




v=0 J=27-26 


343477.843 


343478.4(10) 


0.67(5) 


15.2(8) 
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Table 5. KC1 line parameters in IRC +10216 



Transition 


v rcst 


y obs 


tt^, — ~, 

jTJdv 


v exp Notes 








(MHz) 


(MHz) 


Kkms-' 


km s" 1 










K 35 C1 






v 


=0 1= 


11-10 


84562.590 


84562.6(8) 


0.14(3) 


13.5(10) 


V 


=0 J= 


12-11 


92246.501 


92246.6(6) 


0.20(2) 


14.5 || 


v 


=0 J= 


= 13-12 


99929.473 


99928.8(10) 


0.21(3) 


14.0(8) 


V 


=0 J= 


=14-13 


107611.429 


107612.6(10) 


0.32(5) 


12.2(10) $ 


V 


=0 1= 


17-16 


130650.412 


130650.3(8) 


0.35(6) 


14.2(8) 


V 


=0 J= 


19-18 


146003.215 


146001.9(10) 


0.37(6) 


11.3(10) 


V 


=0 1= 


20-19 


153677.427 


153677.8(15) 


0.34(8) 


15.5(15) ! 


V 


=0 J= 


23-22 


176690.360 


176692.2(15) 


0.62(10) 


16.0(15) 


v 


=0 J= 


27-26 


207348.712 


207347.8(10) 


0.48(5) 


11.3(10) 


V 


=0 J= 


28-27 


215008.213 


215008.8(6) 


0.40(3) 


14.7(5) 


V 


=0 1= 


29-28 


222665.524 


222665.3(5) 


0.68(3) 


12.1(6) 


V 


=0 1= 


30-29 


230320.564 


230320.5(8) 


0.60(3) 


13.2(6) 


v 


=0 J= 


31-30 


237973.257 


237972.0(10) 


0.46(4) 


13.6(8) 


V 


=0 1= 


-33-32 


253271.284 


253271.1(5) 


0.90(4) 


14.0(6) 


V 


=0 1= 


=34-33 


260916.462 


260916.2(10) 


0.58(4) 


11.8(8) 


V 


=0 1= 


=36-35 


276198.754 


276198.7(15) 


0.54(6) 


12.3(15) 


V 


=0 1= 


=37-36 


283835.711 


283838.0(20) 


0.80(3) 


13.2(5) 


V 


=0 J= 


39-38 


299100.855 


299105.6(20) 


0.71(8) 


16.7(20) 


V 


=0 1= 


40-39 


306728.885 


306727.9(10) 


0.85(7) 


15.7(10) 


V 


=0 1= 


41-40 


314353.781 


314354.4(10) 


0.66(5) 


15.3(10) 


V 


=0 1= 


=42-41 


321975.467 


321975.4(15) 


0.34(5) 


7.1(15) 


V 


=0 J= 


=43-42 


329593.862 


329594.5(10) 


0.62(6) 


10.5(10) 


V 


=0 1= 


=44-43 


337208.889 


337209.4(10) 


0.40(3) 


13.6(10) 


V 


=0 1= 


45-44 


344820.468 


344818.4(20) 


0.66(9) 


17.5(20) 










K J/ C1 






V 


=0 1= 


-11-10 


82159.153 


82159.0(10) 


0.08(2) 


14.5 || ! 


V 


=0 J= 


= 12-11 


89624.771 


89625.0(10) 


0.08(2) 


14.5 || S ! 


V 


=0 1= 


= 13-12 


97089.504 


97089.0(7) 


0.08(2) 


13.3(8) 


V 


=0 1= 


= 14-13 


104553.276 


104553.3(15) 


0.11(3) 


14.5 || $ ! 


V 


=0 1= 


= 15-14 


112016.016 


112017.1(15) 


0.06(2) 


14.5 || $ ! 


V 


=0 1= 


28-27 


208903.957 


208903.2(15) 


0.18(2) 


14.5 || ! 


V 


=0 1= 


30-29 


223782.837 


223782.8(10) 


0.19(2) 


10.4(10) 


V 


=0 1= 


=33-32 


246084.350 


246084.2(10) 


0.12(2) 


11.0(10) ! 


V 


=0 1= 


=35-34 


260939.958 


260942.4(20) 


0.42(8) 


12.0(20) 


V 


=0 1= 


-36-35 


268363.919 


268365.2(20) 


0.27(4) 


14.5 || S 


V 


=0 1= 


39-38 


290619.543 


290619.2(20) 


0.42(6) 


11.6(20) 


V 


=0 1= 


=40-39 


298032.418 


298034.7(20) 


0.19(3) 


7.1(15) 


V 


=0 1= 


=43-42 


320253.005 


320251.0(20) 


0.22(3) 


12.0(20) 


V 


=0 1= 


45-44 


335050.955 


335051.0(20) 


0.14(2) 


6.4(20) ! 
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Table 6. A1C1 line parameters in IRC +10216 



Transition 


v rcst 


y obs 


tt^, — ~, 

jT'rfv 


v exp Notes 




(MHz) 


(MHz) 


Kkms-' 


km s" 1 






A1 35 C1 






v=fj J=6-5 


87458.226 


87458.1(3) 


0.53(4) 


14.5(8) $ 


v=0 J=7-6 


102031.870 


102031.6(2) 


0.90(4) 


14.9(6) 


v=0J=10-9 


145744.527 


145744.6(5) 


2.85(10) 


14.2(7) $ 


v=0 1=11-10 


160312.055 


160312.0(2) 


3.31(7) 


14.0(4) 


v=0 1=12-11 


174877.605 


174877.4(4) 


2.21(10) 


13.9(8) $ 


v=0 J=14-13 


204002.051 


204002.0(4) 


3.13(5) 


14.0(6) 


v=0 1=15-14 


218560.586 


218560.5(5) 


3.20(6) 


13.8(5) 


v=0 J=16-15 


233116.423 


233115.7(5) 


5.2(1) 


14.2(4) 


v=0 1=17-16 


247669.382 


247669.4(5) 


6.8(2) 


14.2(5) $ 


v=0 1=18-17 


262219.283 


262219.0(8) 


6.6(2) 


14.8(8) $ 


v=0 1=19-18 


276765.946 


276767.1(6) 


5.9(1) 


14.1(5) 


v=0 1=20-19 


291309.191 


291309.7(6) 


6.6(1) 


14.3(5) 


v=0 1=21-20 


305848.837 


305848.9(5) 


6.4(1) 


14.4(5) 


v=0 J=22-21 


320384.706 


320384.1(6) 


5.8(1) 


13.2(7) 


v=0 1=23-22 


334916.617 


334916.7(5) 


4.45(10) 


13.7(6) 


v=0 J=24-23 


349444.289 


349443.3(10) 


7.4(2) 


14.5 || $ 






AP'Cl 






v=0 J=6-5 


85403.976 


85403.4(8) 


0.11(2) 


14.5 || 


v=0 J=7-6 


99635.372 


99635.5(4) 


0.35(3) 


14.6(7) 


v=0 J=8-7 


113865.568 


113865.4(7) 


0.29(3) 


13.3(10) 


v=0 J=10-9 


142321.673 


142322.5(10) 


0.84(8) 


10.9(10) 


v=0 1=11-10 


156547.240 


156547.1(10) 


1.77(15) 


14.5(15) $ 


v=0 1=15-14 


213428.929 


213428.6(4) 


2.19(5) 


14.0(5) 


v=0 1=16-15 


227643.350 


227643.3(4) 


2.40(4) 


14.3(4) 


v=0 1=17-16 


241855.027 


241855.2(3) 


2.26(4) 


13.8(3) 


v=0 1=18-17 


256063.788 


256063.8(4) 


1.52(3) 


13.8(4) 


v=0 J=19-18 


270269.462 


270269.2(5) 


2.53(3) 


14.2(5) 


v=0 1=20-19 


284471.877 


284473.1(8) 


2.66(3) 


14.0(5) 


v=0 J=21-20 


298670.862 


298672.0(8) 


1.79(2) 


13.9(5) 


v=0 J=22-21 


312866.245 


312865.3(10) 


1.79(3) 


14.5 || S 


v=0 1=23-22 


327057.854 


327057.1(8) 


2.46(4) 


13.8(8) 


v=0 1=24-23 


341245.517 


341245.2(6) 


1.66(3) 


13.2(6) 


v=0 J=25-24 


355429.063 


355430.9(10) 


1.68(4) 


14.2(8) 



Table 7. A1F line parameters in IRC +10216 



Transition 


y rest 






Vexp 


Notes 




(MHz) 


(MHz) 


Kkms- 1 


km s -1 








A1F 








v=0 J=3-2 


98926.768 


98926.8(5) 


0.93(8) 


14.1(7) 


s 


v=0 J=4-3 


131898.841 


131899.1(15) 


1.95(15) 


17.0(15) 




v=0 J=5-4 


164867.899 


164867.4(8) 


2.8(1) 


12.8(8) 




v=0 J=6-5 


197833.190 


197833.1(3) 


2.6(1) 


14.4(4) 




v=0 J=7-6 


230793.958 


230793.6(3) 


3.5(1) 


14.3(5) 




v=0 J=8-7 


263749.452 


263749.3(4) 


5.2(1) 


14.5(4) 




v=0 J=9-8 


296698.916 


296699.0(3) 


4.7(1) 


13.8(4) 




v=0J=10-9 


329641.598 


329640.0(8) 


5.1(1) 


14.3(5) 
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Table 8. NaCN line parameters in IRC +10216 



Transition 


^rest 




J L A av 


v exp 


Notes 




(A/HI '7 A 


^ivinz ) 


JV Kill S 


km s 






in av_iN 










S A 

3l,4-*l,3 


808AA Q01 


OUO't0.7^j ) 


70^ i 
U. / y\D ) 


L j. lyj) 




A S 


QO^QA 




o Q■^^'S^ 
u.yj^j ) 


1 A Sl(/A\ 
14.0^4 ) 




A S 

"0,6-30,5 


CH70A OQ7 
y jzuo.uyz 


Q^70A (\{^\ 


1 1 0^7^ 


1 A KA\ 
LH.Z.yH ) 




A S i a z 

°5,2-->5,l + 05,1-35,0 


Q^A^7 ^83 _i_ Q3A^17 ^8A 
7j0j / .joj t yjOj /.Jo4 


Q^A^7 1^10^ 


U. 14^Z ) 


1 A AH OA 
14.4^1U^ 


£■ 


A S 

02,5 "3 2,4 


0A71 7 ^A^ 
7J / 1Z. j'+j 


Q^7 17 AH (Y\ 
yj 1 LZ-.t+y LVJ) 


o ash ^'i 

U.Dj^ 1 J ) 


1 A ^ 
14. J 


II *t 
II M* 


A S _i_ A S 

D4,3-J4,2 + 4,2-34,l 


Q^7^8 Q70 _i_ Q^7QQ 07^ 

7j / jo.? /u t yj / jy.uzj 


yj 1 jo. Dy^+j 


O ^^(A\ 
U.J J ^4 ) 


1 7 7^8^ 
1Z. iyo) 


i? 
X, 


A S 

03,4-33,3 


yjo j 


7JOJO.J^ / ) 


AA^'A^ 


1 -2 o//;\ 
1 j.o^O ) 




A S 

03,3-33,2 


Q^8AS A7^ 
yjo4o.OZj 


Q^SzlS 8^7^ 
7jo4o.oV / ) 


O AA(&\ 
U. 44(^0^ 


1 A ^/AA 
14. J^O ) 




A S 

02,4-32,3 


y4jj4.oUj 


Q/1QQ/1 Q/'Q\ 

V4jj4.y(j ) 


1 (Y)(f^\ 
l.Uz(D) 


1 /I A//1A 
14.0(4) 




A S 

"1,5-31,4 


QAQ^Q SOS 

yoyjy.ouo 


QAQ^Q Qf1\ 
yOyJy.yyj ) 


L .UU(,D ) 


1 A O/'AA 
14.U^4 J 




7 A 
'1,7-01,6 


1 O^Q^ 77^ 
1UJ jy j.Z 1 J 


1Ujj7J.U(,j ) 


1 7A('8 1 l 
1 .Z4^o ) 


1 A 7/^A 
14. /pj 




7 A 

'0,7-Oo,6 


108A71 QS7 
!Uo4 / 1 .70 / 


1 OSA77 0^7^ 


1 77l'7^ 
1 .Z /(^ 


1 A 7/^ 
14. /^4J 




7 A _i_ 7 A 

'5,3-05,2 + '5,2-05,1 


100,7^0 1^7 _i_ 1007^0 1^3 

lUyZjU.ljZ + lUyZjU.ljj 


1 OQ7AQ f*(£*\ 
LVyZ^+y.OyO ) 


O ^lAf^ 
U. j4^Z ^ 


1 A 7^8^ 
14.Z^o ) 


X. 


7 A 

'2,6-02,5 


10Q781 A78 


10Q781 ^fSTl 


U.oD^o ^ 


1 A O/'AA 




7 A -i-7 A 
/4,4-04,3 + '4,3-04,2 


10Q^7^ 77A _i_ 100A7^ A^8 
lUyj/j.Z/0 t IUVj /D.4JO 


1 OQ^7^ 7^7^ 


O AAf'A 1 ! 


i j.zp J 


£ 
X, 


7 A 

'3,5-03,4 


100^01 qao 
iuvjui .you 


100^01 Q/'A^ 
LVyJVL .y{0 ) 


U.OZ^J ) 


1 A ^ 
14. j 


II <t 
II 


7 A 

'3,4-03,3 


1 OCK7A 8A3 

i uy j Z4. ot'J 


1 OQ^7^ ^/"A^ 


O ^71^ 


1 A ^ 
14. j 


II t 
II * 


7 A 

'2,5-02,4 


1 1 07AQ Ql 7 


1 1 0770 f\f1\ 
L 1UZ /U.U^j ^ 


i 1Q/'A^ 


1 A Q/'^A 
14. yyj ) 


>> 


7 A 
'l,6-Ol,5 


I 1 ^OAO 7 1 Q 

I I JU4U.Z1;? 


1 1 ^ozin 1 /"A^ 




1 A 'J/'AA 
14. j^O J 


>> 


Q 7 
01,7- '1,6 


17Q081 "307 

izyuo 1 . JUZ 


17QOQ1 A/'QA 


1 77/1 f\\ 
L . / Z{ 1 U ) 


i q en OA 
1 J.O^IU; 




Q S 

71,9-01,8 




1 Q^07 A(L(X\ 


1 77/1 7A 
1 .ZZ^ LI) 


i j.y^i j) 




Q 8 

"0,9-O0,8 


1 ^sa^7 ock 
1 joojz.uyj 


1 ^8A^7 {^(^ 


1 ^i'A'* 
1 . J j(,D^ 


i j.yp j 




74,6-04,5 + 7X5-64,4 


1 AOAAA AOS _i_ 1 AOAA7 807 


1 AOAAA 'Xf 1 7^ 


O A7('8'l 
U.D f\o) 


1 a ah OA 

1 J.411U,) 


£ 


Q 8 

"2,7-02,6 


1 A7zll fl S77 


1 ZL7A1 zl/'S^ 


SAi"^ i 
u.oop ^ 


1 A ^/"AA 
14. jyH ) 




Q 8 

71,8-01,7 


14DU / J. JOD 


1/1^07^ ^/"^^ 


1 77('A^ 


1 ^ S/'AA 
1 j.O^O ^ 


4 


1 n q 

fJo.10-70,9 


1 ^^^^7 Q77 


1 ^^^^7 f\( 1 O^ 
Ijjjj /.U^IUJ 


1 l^i^ 


ic 7/1 OA 
1 J.Z^IU; 




1 n q 

102,9-72,8 


1 ^8^7 Q7^ 


L J jo jo. 1 ) 


i OA^'a^ 

1 .U0(,o ) 


1 ^ An OA 
1 J.4^1U; 




1 n q x in q 

104,7-^4,6 + iU4,6-74,5 


1 ^A^7^ _l 1 ^^^7^ Q^n 
IjO^Zj.jj'4 + LJDjZj.yjV 


1 ^A^7zl ^(\ (Yi 
1 j0jZ4.j(1U,J 


1 AA/" 1 O^ 
1 .40^ LU) 


1 a 7n OA 

1 J. / ^1U; 


£ 
X. 


1^3,8-73,7 


1jDj £ t1 1 JjU 


1 jOj4U.7^1 J ) 


AAfQ i 


1 A OA 
14. J^IU^ 


*t 


1 n o 

103,7-73,6 


1 jOOo4.jjo 


1 J00o4. i\o) 


O SAf'A^ 
U.oO^O ) 


I 1 0/"AA 

I I .U^O J 




1^2,8-72,7 


1 ^SA 1 A 7A7 


1 ^8A1 7 7/'8^ 


1 07('8 1 i 
1 .U i\o) 


1 A S/'SA 
14.0^0 J 




in o 

1( Jl,9-71,8 


1 A1 H 1 A 1QA 
IDIUI^-. ly^ 


1 A 1 01 zl S/'^A 
101U14.o^J ^ 


1 . 1 J^J ) 


1 A 1 /^A 
14.1pJ 




11 in 1 1 1 in 

115,7-11)5,6 + 115,6-105,5 


1717^^ ASIA -j_ 1717^^ ^AC\ 


1 7 1 7^ Ai' 1 f\\ 
L 1 L 1 jj.4^1UJ 


o son 7^ 

V.oyyLZ) 


Q 1 n 7A 
y.L{LZ) 


x- 
X. 


11 in 

112,9-1U2,8 


1 7zlQHA 1 7^ 


1 7zlQ0^ 

L 1 ^yxjj.DyJ ) 


U.4o^j ^ 


1 ^ ^fAA 
1 J . J^4 J 




11 in 

111,10-101,9 


1 /Oooy.yyo 


1 7A88Q 1 ^8^ 
1 lOooy.Lyp) 


O AAf'A'l 


1 ^ A/'SA 




19 11 

1^1,12-1 ll.ll 


17QQ71 QOA 
1 lyyLY .y\JO 


17QQ70 8M7^ 
1 lyyZXf.oyLZ) 


QAH Ci\ 
\J.y\j[ L\J) 


1 7 AH SA 
IZ.O^l J 7 




1 ^ 19 

1 30,13-14),12 


1 Q7M 1 AA7 


1 Q7A1 1 S/'A^ 
Ly /Dl 1 .0^0 ^ 


O i ^'?^'^^ 
U.jjp ) 


1 ^ 0/"7A 
1 J.U^ / ) 




1 1 19 

-132,12-1^2,11 


707 1 fH 717 
ZUZ1UJ. / 1 / 


707 1 OA 0^70^ 

zuziuo.y^zu^ 


1 OOf^ 1 ! 
1 .UV^J ) 


1 A OOOA 
lO.U^ZU^ 




1 ^ 19 1 11 19 

136,8-1A,7 + 136,7-1A,6 


ZUZDoo.Joj + ZUZO00.J0/ 


707A8^ 0^70^ 
ZUZOoJ.y^ZUJ 


U.0j(^j ) 


1 a on ^A 

14.U11 J ; 


£ 
X. 


1 ^ 19 -1- 1 ^ 19 

135,9-1^5,8 + 135,8-1/5,7 


70^000 52^8 -1- 70^001 1 ^1 
ZUjvJUU.ojo t ZUjUUI.IjI 


70^007 8^1 Ci\ 
ZUjUUZ.O^lU^ 


Ql fAi 

u.y LyO) 


1 o on sa 

1 \J.yy 1 J ^ 


J? 

X. 


1 ^ 19 

133,11-1^-3,10 


70^^77 1 ^7 
ZUODZ / . 1DZ 


70^A77 8^1 f\\ 
ZUjOZ / .ol^lU^ 


U.oZ^O^ 


ic 7/1 OA 
1 J.Z^IU; 




1 ^ 19 

133,10-1^3,9 


7Hzl 1 7^8 
ZU41 JO. / JO 


70A 1^8 Ai' 1 O^ 
ZU41 JO.D^IU^ 


O A s s(A\ 
U.4j^4 ) 


11 An OA 
1 1 .O^IU; 




1 ^ 19 

132,11-1^2,10 


7H7A^9 7^7 
ZU /Ojo.Zj / 


707A^7 (\( 1 f\\ 
ZU /Oj / .U^IUJ 


U.oZ^OJ 


1 a An OA 
1 j.O^IU; 




1 ^ 19 

131,12-1^1,11 


ZUo4Uo.UU7 


708A08 A/'AA 


O 7A('^^ 
U. /0(,J ) 


1 ^ ^AA 
1 j. J^O ^ 




1 a n 

1*1,14-131,13 


70QAQ^ 


70QAQ7 A^l 
ZUV4y / .4(1 J ) 


U. J7(,J ^ 


1 1 *;n ^a 
1 1 . jii j ) 




1 A 1 ^ 

1*0,14-130,13 


71717/I 07Q 
Z 1Z IZH-.UZ? 


71717^ t ( t 0^ 
ZlZlZj. L\L\j) 


u.jzp ) 


1 ^ An OA 
1 J.4(1U,J 




1 A 1 ^ 

1*3,12-133,11 


7 10^71 OA8 
Ziyj>Zl .U4o 


7 1 Q^7^ *,(\ C\\ 
LLyoLD. j\L\j) 


O AA(Si\ 
U. 44(^6 ^ 


1 a 1 n OA 
1 j. I^IU; 




1 A 1 ^ 
1+3,11-133,10 


77008A 


77008A 7^1 
ZZUUoO.Z^l J ) 


O ^7/ A^ 
U. JZ^O^ 


1 ^ 7^1 OA 
1 J.Z^IU; 




1 A 1 ^ 
1*1,13-131,12 


77AO^O A^8 


77AO^O 
ZZ4UJU. J{J ) 


O '^7^'A^ 
U. JZ^j 


1 7 A/"AA 
1Z.4^0 ) 




1 A 1 ^ 

1*2,12-132,11 


77/1087 Al A 
ZZH-UoZ.D 1 D 


77/1087 0{ c \\ 
ZZH-UoZ.Z^J ) 


A7('^ i 
U.OZ^J ) 


1 ^ 7/'AA 
1 j.Z^O ) 




1 S 1 A 

131,15-1*1,14 


77zl7^7 7^7 
ZZ4ZJ3Z.ZJ / 


77A7^7 1 (f*\ 
ZZ4ZjZ. 1^0^ 


O ^ 1 ^'A^ 
U.J Ly*) 


1 A 1 /"AA 
14. 1^0 J 




1 S 1 A 

130,15-1*0,14 


77A^81 1^7 
ZZO JO l.lj/ 


77A^SO 'xn 0^ 
ZZOjoU. j^IUJ 


A7^'A^ 
U.D / ^4^ 


1 ^ ^fAA 




1 S 1 A 

132,14-1*2,13 


7^77A8 1 7Q 
ZjZ / 4o. IZ7 


7^77A7 Ai' 1 f\\ 
ZjZ / 4 / .O^IU^ 


O A8^'A^ 


1 A 1 /"SA 
14. 1^0 J 




IS 1 A , 1C 1 A 

137,9-1*7,8 + 137,8-1*7,7 


7^1^ 1 O ^8^ _i_ 7^1^ 1 n ^8^ 


7^^ 1 o on C\\ 
ZjjjIU.U^IUJ 


O ^Oi^ 
U. jU(,j ) 


1 J.U^IU; 


£ 
X. 


1 S 1 A ■ 1 S 1 A 

136,10-1*6,9 + 136,9-1*6,8 


ZJJ07J.JUU + ZjJ07J.jlj 


7^^8QA 7n n\ 

Zjjoyo.z^iu^ 


U. Dyyj ) 


1 o on OA 
lu.y^iu; 


£ 


1 S 1 A _i_ 1 S 1 A 

135,11-1*5,10 + 135,10-1*5,9 


7^zl78Q I'XA _i_ 7^1/1700 771 
Zj4Zo!7. / J4 + ZJ4ZVU. / / 1 


7^A7QO A( 1 C\\ 


O ^Of^ 


1 7 A/'SA 
IZ.O^o ) 


£ 
X. 


1 S 1 A 

134,12-1*4,11 


Zj4 /jo.1 oj 


Zj4 / j / . j^o ^ 


U. jU(,j ) 


1 7 Q/'SA 




1 ,-,-14, , , 

^-m,^ ^-^3,11 




236079 6(^ 


0.20(2) 


11.1(10) 




15 U 4-14 U3 


239546.409 


239546.2(6) 


0.58(4) 


12.5(6) 




15 2 ,i3-14 2 ,i2 


240507.236 


240506.8(5) 


0.58(4) 


13.4(6) 




162,15-152,14 


248003.167 


248002.7(5) 


0.55(4) 


13.3(6) 




167,10-157,9 + 167,9-157,8 


249079.914 + 249079.915 


249079.1(8) 


0.49(6) 


14.2(8) 


£ 


166,11-156,10 + 166,10-156,9 


249502.758 + 249502.785 


249503.2(8) 


0.64(7) 


12.9(8) 


£ 


165,12-155,11 + 165,11-155,10 


249943.220 + 249945.104 


249944.0(6) 


0.69(7) 


12.4(6) 


£ 


17 U7 -16 U 6 


253611.361 


253610.6(8) 


0.54(6) 


12.8(8) 
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Table 8. Continued. 



Transition 






v obs 


i'T'dv 

J A 


Vexp 


Notes 






(MHz) 


(MHz) 


K km s~ l 


km s _1 




1 ' U,l / iV/ U,lo 




255396.376 


255397.6(8) 


0.78(7) 


14.0(5) 




1 / ih t- nin 4- 1 / i f\ o 
A ' 1U,/ lv -'lU,o ' 1 ' 1U,0 




263154.449 + 263154.449 


263155.7(20) 


0.36(5) 


18.4(20) 


£ S 


' ' 2,iO iv/ 2,15 




263210.770 


263210.4(8) 


0.38(4) 


14.0(8) 




1 /n n- f"ln o t™ 1 /no — 


1607 

±v/ y, / 


263682.103 + 263682.103 


263683.6(10) 


0.18(2) 


12.4(10) 


£ 


1 / o in~ 1 t)o ^ ~r~ 1 / o n — 

1 ' o,lU A yj o,y ' 1 1 x,y 


1 (So Q 


264178.070 + 264178.070 


264177.2(10) 


0.41(3) 


15.8(10) 


£ 


1 If. 11-1 Oft 11 + l /ft 1 1 


" 1 ^6, 10 


265112.789 + 265112.843 


265111.1(15) 


0.66(4) 


15.3(15) 


£ 


175^3-16512 1^5,12 


■16s 1, 


265603.232 + 265606.523 


265606.4(15) 


0.74(5) 


13.2(15) 


£ 


1^4,14" 1"4,13 




266168.587 


266167.4(15) 


0.46(5) 


17.8(15) 




1 1 A IT 1 tl^ IT 

1 '4,13 AV, 4, 12 




266283.597 


266281.4(20) 


0.39(6) 


14.5 


II $ 


1 fi 1 c- 1 riM i ^ 

A ' 3,15 AV7 3,14 




266346.964 


266347.0(10) 


0.61(5) 


14.1(10) 




1 8 1,18" 1^1,17 




268257.278 


268259.9(20) 


0.46(5) 


13.4(15) 


$ 


1^3,14" ^^3,13 




268287.549 


268288.5(10) 


0.36(4) 


12.3(10) 




1 On i o — 1 1 (\ \n 

1 U U, Ifi - 1 - ' U, 1 / 




269780.806 


269780.0(8) 


0.48(4) 


11 7(8) 




1/i ift-1 Oi n 
1 ' l,lo AV, 1,15 




270217.365 


270218.8(10) 


0.58(4) 


12.7(8) 




1 /IK- fll 1 ,1 

A ' 2, 15 A "2,14 




273268.589 


273270.7(15) 


0.50(6) 


10.1(15) 




1 O i h — 1 / o i a 

1 "2,1 / A ' 2,16 




278369.353 


278368.5(10) 


0.42(4) 


15.1(10) 




1 on i n — 1 / n n ~r~ 1 on n — 

A "y,iu A ' y,y 1 A "y,y 


179 8 

A ' y,s 


279177.286 + 279177.286 


279173.8(20) 


0.26(4) 


10.5(20) 


£ 


loo ii — l / o in ~r~ loo in 

L "0,1 1 1 ' 8,1U 1 AVJ o,lL 


-17l!9 

' ' s,y 


279710.006 + 279710.006 


279715.1(20) 


0.44(6) 


14.5 


£ II $ 


1 0*7 1 T — 1 /Til 4- 1 m 1 1 

A "/,12 ± ' /,11 T A "/,ll 


"1^7,10 


280219.750 + 280219.751 


280220.6(10) 


0.49(5) 


14.6(10) 


£ 


1 $4,14" 1^4,13 




282067.503 


282070.2(20) 


0.43(5) 


15.2(15) 




1 9[ 5 19 - 1 8 1,18 




282876.644 


282878 9f20) 


0.47(5) 


13 5(15) 




1 7n in— 1 on 10 




284161.712 


284161.3(10) 


0.43(3) 


11 4(8) 




1 fS 1 1 C — 1 / 1 I A 

± "3,15 A ' 3,14 




284514.720 


284516.3(15) 


0.38(4) 


10.8(10) 




1 rS i i h — 1 / i i a 

A "1,1 / A ' l,lo 




285355.280 


285353.0(20) 


0.45(5) 


15.3(15) 




1 iSo i^—l / 1 i £ 

- 1 "2, lo A ' 2,15 




289567.202 


289566.9(10) 


0.48(4) 


13.8(10) 




1 ^2, 18" 1 ^2,17 




293477.704 


293481.9(20) 


0.77(8) 


16.7(15) 




1 ^n i i — 1 rSn i n 4~ 1 ^n i n 

A - y y,n iu y,iu ^ ^-'y,!!. 


18 

- 1 09,9 


294669.796 + 294669.796 


294669.1(20) 


0.33(5) 


13.2(15) 


£ 


1Q 7 ,,-18t n + 19t n 


" 1 "7, 1 1 


295790.308 + 295790.311 


295788.6(20) 


0.30(5) 


14.5 


£ II $ 

II 4" 


19, 1,-185 ii + 19s ii 

a>, 5,1j A "5,14 ^ A -0,14 


" 1 05, 13 


296944.126 + 296953.235 


296948.3(20) 


0.28(4) 


15.0(20) 


£ 


20, 9n-19i iq 

^"i,2u A -^i,iy 




297471.433 


297470.7(20) 


0.24(4) 


11.0(20) 




Aj/ 4,15 A "4,14 




297879.787 


297879.7(15) 


0.20(3) 


12.0(15) 


$ 


20o,20~ l^ojy 




298544.239 


298542.5(20) 


0.28(4) 


11.2(15) 




1^1 1 O — 1 O 1 IT 

i -'l,lo 1(J 1,1/ 




300354.316 


300354.0(10) 


0.44(4) 


11.7(10) 




■^■"o, 13 1 - , o,12 1 ^-"8,12 


■19s ,, 


310769.695 + 310769.695 


310767.2(20) 


0.22(3) 


12.6(20) 


£ 


20t 1/1-19? 1 2 + 20t n 

^-"/,14 J -- y /,13 ^ ^-"/,13 


■19v ,? 

J --' /,12 


311361.332 + 311361.337 


311361.8(10) 


0.20(2) 


14.2(10) 


£ 


20fi if-19fiid + 20^ ia 

"o, 1 j ij, o,14 1 ^-"0,14 


■19fi 13 

1 -'D,13 


311961.365 + 311961.708 


31 1961.6(10) 


0.35(3) 


15.0(10) 


£ 


/I i t i — / 1 1 1 on 

A 1,21 ^-"1,20 




312043.593 


312043.2(8) 


0.33(2) 


12.0(8) 




2, 1 n t i — Z^l /n on 
^ A U,21 "11,21) 




312931.171 


312930.1(15) 


0.13(1) 


12.6(15) 




20i ,s-19i n 




313196.095 


313198.4(20) 


0.12(2) 


14.5 


II 
II 


20i n-19zi i,; 




313369.371 


313371.4(15) 


0.22(2) 


16.6(15) 




204 ,fi-19d is 




313724.992 


313721.0(20) 


0.44(5) 


14.5 


|| 
II 


20i ,o-19, ,» 

z - vy i, iy A -'1,1s 




315213.834 


315212.9(15) 


0.39(4) 


13.9(10) 




20i n-19i if, 




317220.575 


317219.5(10) 


0.38(3) 


12.8(10) 




2 1 1 20 - 20i 19 




329937.919 


329938.4(20) 


0.22(3) 


11.4(20) 




23l,23-22i,22 




341127.341 


341128.1(15) 


0.37(3) 


14.8(15) 




23o,23-22 ,22 




341721.282 


341717.9(20) 


0.30(4) 


12.3(20) 




228,15-218,14 + 22s,14 


"218,13 


341823.295 + 341823.295 


341823.2(10) 


0.40(4) 


12.7(10) 


£ 


226.17-210JO + 22 6 j 6 


■21 6 ,15 


343211.017 + 343212.019 


343213.0(20) 


0.42(4) 


15.2(20) 


£ 


223,20-213,19 




344269.527 


344267.6(20) 


0.19(4) 


11.5(20) 




22i,2i-21i,2o 




344535.644 


344539.8(20) 


0.22(3) 


10.4(20) 




224,19-214,18 




344854.243 


344857.9(20) 


0.25(4) 


10.2(20) 





